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ABSTRACT 



Various elfective temperature scales have been proposed over the years. Despite much work and the high internal precision usually 
achieved, systematic differences of order 100 K (or more) among various scales are still present. We present an investigation based 
on the Infrared Flux Method aimed at assessing the source of such discrepancies and pin down their origin. We break the impasse 
among different scales by using a large set of solar twins, stars which are spectroscopically and photometrically identical to the 
Sun, to set the absolute zero point of the effective temperature scale to within few degrees. Our newly calibrated, accurate and precise 
temperature scale applies to dwarfs and subgiants, from super-solar metallicities to the most metal-poor stars currently known. At solar 
metallicities our results validate spectroscopic effective temperature scales, whereas for [Fe/H] < -2.5 our temperatures are roughly 
100 K hotter than those determined from model fits to the Balmer lines and 200 K hotter than those obtained from the excitation 
equilibrium of Fe lines. Empirical bolometric corrections and useful relations linking photometric indices to effective temperatures 
and angular diameters have been derived. Our results take full advantage of the high accuracy reached in absolute calibration in recent 
years and are further validated by interferometric angular diameters and space based spectrophotometry over a wide range of effective 
temperatures and metallicities. 
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1. Introduction 

The determination of effective temperatures (reff) in F, G and K 
type stars has a long and notable history. Because of their long 
lifetimes these stars retain in their atmospheres a fossil record of 
the chemical elements in the interstellar medium at the time of 
their formation. The stellar effective temperature is of paramount 
importance for reliable abundance analyses and thus for improv- 
ing our understanding of Galactic chemical evolution. 

Stellar abundances are now routinely derived from high res- 
olution spectra, model atmospheres, and spectrum synthesis. 
, While each of these ingredients have their own issues regard- 
■ ing systematic uncertainties, the dominant source of error is in 
' many cases the adopted T^ff of the star Several indirect meth- 
ods of Teff determination have been devised to avoid the com- 
plications introduced by the measurement of stellar angular di- 
, ameters, w hich are necessary to derive T^m from basic prin- 
ciples (e.g. iHanburv Brown et al]|1974t Ivan Belle & von BraunI 
l2009l) . Thus, most published values of Teff are model-dependent 
or based on empirical calibrations that are not free from system- 
atics themselves. 

It is therefore not surprising to find discrepancies among 
published Teff values. The ionization and excitation balance of 
iron lines in a ID LTE analysis is routinely used to derive effec- 
tive temperatures as well as log^ and [Fe/H]. While for a sam- 
ple of stars with similar properties this method can yield highly 
precise r elative physical parameters (Melen dez et al.l l2009at 
iRamfrez e t al. 2009, see Section [3] for its use on solar twins), 
non-LTE effects and departures from homogeneity can seriously 
undermine effective temperature determinations, especially in 



metal-poor stars (e.g. 'Aspl undll2005h . Similarly, the line-depth 
ratio technique has high internal precision, claiming to res olve 
temperature differences of order 10 K ( e.g. iGrav & JohansonI 
119911: lGravlll994l: iKovtyukh et aP l2003h b i it it is not entirely 
model independent (e.g. lCaccin et alj|2002l : iBiazzo et"al]|2007l) 
and the uncertainty on its zero point can be considerably large. 
Another popular method for deriving T^ff in late-type stars is 
provided by the stud y of the hydrogen Balmer lines, in p artic- 
ular Ha and H/3 (e.g. iNissen et al.l l2007l 'Fuhr mannll2008 ). For 
H lines uncertainties related to o bserv ations and line broadening 
(Barklem et al. 2002), non-LTE ("Barklem 2007) and granulation 
effects (Asplund 2005; Ludwig et al. 2009; Per eira et al. 20091) 
all influence the estimation of effective temperatures. 

In such a scenario, an almost model independent and ele- 
gant technique for determining effective temperatures was in- 
troduced in the late 70's by D. E. BlackweU and collab orators 
dBlackwell & Shalhsl [19771; iBlackweU et all 1 19791 [T980h under 
the name of InfraRed Flux Method (hereafter IRFM). Since 
then, a number of authors have applied the IRFM to deter- 
mine effective temperat ures in sta rs with different spectral types 
and rnetalhcities (e.g. 



Bell & G ustafsson 
rrr 



1989; Alonsoetal 



'1996a; Ramirez & Melendez 2005a; Casagra nde et ah i200 
Gonzalez Hernandez & Bonifacio 2009). The main ingredient of 
the IRFM is infrared photometry, with the homogeneous and 
all-sky coverage provided by 2MASS being the de facto choice 
nowadays. As such, the IRFM can now be readily applied to 
many stars, making it ideal to determine colour-temperature- 
metallicity relations spanning a wide range of parameters. The 
effective temperatures determined via IRFM are often regarded 
as a standard benchmark for other techniques. Whilst they have 
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high internal accura cy and are essentially free from non-LTE and 
granulation effects ("Asplun d & Garcia Perejl200ll: ICasagrandd 
[2009, Ramirez et al. in prep.), the reddening and absolute flux 
calibration adopted in such a tec hnique can easily introduce a 
systematic error as large as 100 K (ICasagrande et al.li2006.) . 

The effective temperatures of dwarfs and subgiants are still 
heavily debated with various Teff scales behaving very dif- 
ferently depending on colours and metallicities. One of the 
most critical discrepancies occur at the metal-poor end, for 
[Fe/H] < -2.5. In their work on the determination of ef- 
fective temperatures via IRFM, Ramirez & Melendez (2005a) 
found temperatures significantly hotter than those previously 
published, in particular those determined using the excitation 
equilibrium method. Diff'erences up to 500 K for the hottest 
(reff ^ 6500K) most metal-p oor ([Fe/H] < -3.0) stars 
were reported (e.g.. lMelendez & R amirez 20041: iMelendez et al.l 
2006b). In this regime, the recent IR FM investigation by 
Gonzalez Hernandez & Bonifacio! (l2009h still supports a tem- 



perature scale significantly hotter than excitation equilibrium 
and B aimer lines, but ~ 90 K cooler than iRamirez & Melendez! 
(l2005a) . 

The abundance pattern measured in metal-poor stars is im- 
portant for our quest to understand Galactic chemical evolution 
and Big Bang nucleosynthesis: two notable examples are the 
oxygen abundance and the lithium trend with metallicity, both 
of which crucially depend on the adopted Tetf scale. For exam- 
ple, a change of -I-IOOK in T^ff would decrease the [O/Fe] ra- 
tio in turn-off metal-p oor stars by ~0.08 dex w hen using the OI 
triplet and Fell lines (!Melendez et al.! l2006a!) . while the same 
chan ge in Tptf would incre ase the Li abundance by ~ .07 dex 
(e.g. lMelendez & Ramfrej!2004t!Melendez et al.!l20093bh . 

At higher metallicities, which encompass most of the stars 
in the solar neighbourhood, the situation is also uncertain, 
with spectroscopic eff'ective temper atures i n rough agreement 
with the IRFM scale of Casagrand e et al.! (|2006). The latter 
is then about 100 K hotter than the IRFM temperatures of 
!Raimrez & Melendez! (i2005b,) whilst the recent implementa- 
tion of Gonzalez Hernandez & Bonifaci 3 (120091) falls in be- 
tween these two extremes. These difi'erences are somewhat puz- 
zling considering that all recent works on the IRFM have used 
2MASS photometry. Effective temperature calibrations are also 
crucial in the context of deriving reliable colours for theoreti- 
cal stellar models, which ap art from few notable exceptions (e.g. 
!VandenBerg & Clem!l2003r) have to resort entirely to theoretical 
flux Ubraries. 

The aim of this work is to uncover the reason(s) behind such 
a confusing scenario and provide a solution to different IRFM 
effective temperature scales currently available in literature. As 
we discuss throughout the paper, this ambitious task is accom- 
plished by using solar twins which allow us to set the absolute 
zero point of the Teff scale. This result is further validated using 
interferometric angular diameters and space-based spectropho- 
tometry. 

The paper is organized as follows. In Section |2] we com- 
pare the results obtained from different authors, focusing in 
particular on tw o independent im pleme ntations of the IRFM 
teanurez & Melendez!l2005ai: ICasagrande et al.ll2006 ) when the 
same input data are used. This approach allows us to precisely 
identify where different Teff scales originate from. A cure to such 
an impasse is then provided in Section [3] The validation of our 
results, together with the new both precise and accurate effec- 
tive temperature scale are presented in Section|4]to|6] We finally 
conclude in Section|7] 



2. Comparing different versions 

In this paper we use an updated ve r sion th e IRFM implementa- 
tion described in !Casagrande et all (!2006!) to nail down the rea- 
sons behind different Teff scales. Our implementation works in 
the 2MASS system and fully exploits its high internal consis- 
tency thus making it well suited to the purpose of the present 
investigation. The core of t he present study is to carry o ut a 
detailed comparison with the !Ramirez & Melendez! (!2005ah im- 
plementation when the same input data are used. For the sake 
of precision, notice that hereafter, when we refer to a Teif de- 
termined by Ramf rez & Melendez! (12005a) we are referring to 
the effective temperatures determined using that implementation 
and not the original values given in that paper. This is because 
of the updated (and more consistent) input data used here and 
also because some of the stars presented in this work do not 
have IRFM T^ff values published yet. In fact, in order to reveal 
trends with metallicity and/or effective temperature, our sample 
is specifically built to cover as wide a range as possible in those 
parameters (Figure [U. 



2.1. Input sample 

The main ingredient of the IRFM is optical and infrared photom- 
etry. The technique depends very mildly on other stellar param- 
eters, such as metallicity and surface gravity, which are needed 
to interpolate on a grid of model atmospheres (see Section l272l i. 
Below we present the papers from which we gathered [Fe/H] 
and log g for all our stars and we also give references to the pho- 
tometric sources. 



The metal-rich dwarfs come from !Casagrande et al.! (!2006h 
who also provide homogeneous and accurate BV(RI)c pho- 
tometry while additional metal-rich dwarfs and subgiants are 
from .Ramirez & Melendez (2005a). We complement the sam- 
ple with a number of moder ately metal-poor stars from the 
study of 'Fabbian etaP (!2009!) and metal-poor turn off stars 
from Hosford et al.! (!2009!) . To investigate the metal-poor end 
of the Teff scale i n more detail, s t ars with reliable input data 
from .Ramirez & jVIelindez! (!2005a!) . !Bonifacio et alJ (!2007h and 
lAoki et al.! ( !2009!) were added. Finally, to explore for the first 
time the hyper-metal-poor regime via IRFM the subgiants 
HE0233-0343 ([Fe/H] < -4 Garcia Perez et al. 2008) and 
HE13 27-2326 ([Fe/H] < -5 [F rebel et alj 120051: lAoki etll] 
20061 !Frebel et al.!!2008!:!Korn et al..,2009) were included. 

New UBV{RI)c photometric observations for some of the 
metal-poor stars in the aforementioned papers were conducted 
by Shobbrook & Bessell (1999; private communication) and 
are given in Table [1] For the remaining stars, optical Johnson- 
Cousins photometry was taken either from Beers et al. (2007h 
or the General Catalogue of Photometric Data (.Mermilliod et al] 
!I997I) . 

Infrared JHKs photometr y for the entire sample is available 
from the 2MASS catalogue (^Skrutskie et al.! !2006!) which also 
includes the uncertainty for each observed magnitude ("j-", "h_" 
and "kjnsigcom"). The infrared median total photometric error 
of our sample is 0.07 mag. (i.e. "j_"-i-"h_"-i-"k_msigcom"= 0.07) 
and never exceeds 0.14 mag. Such an accuracy in the infrared 
photometry implies a mean (maximum) internal error in Tgff of 
25 K (50 K). Notice that the effective internal accuracy is slightly 
worse because of additional uncertainties stemming from the op- 
tical photometry, [Fe/H] and log^. Altogether our final sample 
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Fig. 1. Distribution of T^g, logg and [Fe/H] for the 423 stars in our sample. 



Table 1. New Johnson-Cousins photoelectric photometry ob- 
tained for some of the metal-poor stars in the sample. Each mea- 
surement comprises an average of 4 observations per star. The 
rms of individual observations are 0.02 for the V magnitude, 
0.015 for the U-B colour and 0.008 mags for B-V, V-R, R-I, V-I 
colours. 



Name U B V Rc Ic 



consists of 423 stars: all have BVJHKs photometry while more 
than half have also (RI)c magnitudes available^. 

Proper reddening corrections are crucial to determine T^ff via 
IRFM. We have tested that 0.01 mag. in E(B - V) translates into 
an IRFM effective temperature roughly 50 K hotter. Reddening 
is usually zero for stars lying \yithin the loc al bubble < 70 pc 
from the Sun (e.g. lLerovlll993l iLallement et all 120031) and so 
we have a dopted E(B - V) - for all stars having Hipparcos 
parallaxes dvan Leeuwenll200"7h and satisfying this requirement 
on the distance. For the remaining stars we updated the redden- 
ing corrections in Ramirez & Melendez (2005a) based on vari- 
ous extinction maps and, in particular for metal-poor stars when 
archive high resolution spectra were available, using interstel- 
lar NaD absorption lines lMelendez et alJ(l2009bl) . In broad-band 
photometry the definition of the effective wavelength of a fil- 
ter (/leff) shifts wi th the colour of the star (e.g. Bessel l et al.l 
1998; Ca sagrande e t al. 2006). Therefore a given E(B-V) colour 
excess must be scaled according to the intrinsic colour of the 
source under investigation. From the reddening E{B - V), we 
computed the extinction in each band adopting the reddening law 
of O'Dpnnell (1994) for the optical and Cardelli et al, (1989) for 
the infrared, using the improved estimation of the stellar intrin- 
sic flux obtained at each iteration to bootstrap the computation 
of the correct /Igff in our IRFM code. 

2.2. The IRFM: pros and cons 

The basic idea of the IRFM is to compare the ratio between the 
bolometric flux fsoiiEsnth) and the infrared monochromatic flux 
?^,^(Earth), both measured at the top of Earth's atmosphere (the 
so-called observational Robs factor) to the ratio between the sur- 
face bolometric flux io-T*^) and the surface infrared monochro- 
matic flux TAigiTss, [Fe/H],logg) determined theoretically for 
any given set of stellar parameters. The latter is called the theo- 
retical Rtheo factor. For stars hotter than about 4200 K, infrared 
photometry longward of ~ 1.2 yum ensures we are working in the 
Rayleigh- Jeans part of a stellar spectral energy distribution, a re- 
gion largely dominated by the continuum which linearly depends 
on Teif and thus only mildly on model atmospheres (Figure 
|2]i. An extension of the technique to cooler effective tempera- 
tures using near-infrared photometry is possible, as shown by 



' Other than being available only for a limited number of stars, we 
did not use U magnitudes because of the little flux emitted in this region 
and the high uncertainties related to the absolute calibration and stan- 
dard ization of this passband in both observed and synthetic photometry 
(e.g. lBessellll200l and references therein). 
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Fig. 2. Top panel: Johnson-Cousins-2MASS filter sets used in 
this work. Middle panel: synthetic solar metallicity spectra at 
different Teff. For the sake of comparison all curves have been 
normalized to unit. Bottom panel: difference in effective tem- 
peratures with - without using {RI)c magnitudes to recover the 
bolometric flux. 



ICasagrande et al.l (l2008l) . but this is outside the purpose of the 
present paper. 

Robs and Rfi„o can be immediately rearranged to determine 
Teff, effectively reducing the entire problem to properly recover 
;FBo/(Earth) and (Earth). Both quantities are determined from 
photometric observations, but an iterative procedure is adopted 
to cope with the mildly model dependent nature of the bolo- 
me tric correction. In our ca se we use the fluxes predicted by 
the lCastelli & Kuruczl (l2004l) grid of model atmospheres starting 
with an initial estimate of the effective temperature and interpo- 
lating at the appropriate [Fe/H] and \ogg until convergence in 
Teff is reached within 1 K. By doing so, we also obtain a syn- 
thetic spectrum tailored to the effective temperature empirically 
determined via IRFM. 

Though we interpolate at the proper [Fe/H] and logg of each 
star, the dependence of the technique on such parameters is mi - 
nor re.g. lRamfrez & MelendeH2005at ICasagrande et al.ll2006h . 
This feature makes the IRFM superior to any spectroscopic 
methods to determine Teff -provided the reddening is accurately 
known- since in the latter the effects of Teff, logg and [Fe/H] 
are usually strongly coupled and the model dependence is much 
more important. 

The errors are estimated using realistic observational un- 
certainties in a Monte Carlo simulation plus the systematics 
arising from the adopte d absolute calibration, as described in 
ICasagrande et alJ (l2006h . With the improved absolute calibra- 
tion used in this paper, systematics amount to 15K in Teff and 
0.3% in bolometric flux (Section [3. 2l i. For stars approximately 
cooler than 5000 K, (RI)c photometry is crucial to properly com- 
pute the bolometric flux. This can be appreciated in the lower 
panel of Figure |2l below this temperature a trend appears us- 
ing BVJHKs magnitudes only. Missing the peak of the energy 
distribution clearly leads one to underestimate the bolometric 
flux thus returning cooler effective temperatures. We have lin- 
early fitted the trend below 5000 K to remove such differences 



in both Teff and J^boI when (RI)c photometry was not available. 
For Teff > 5000 K no obvious trend appears: constant offsets 
of merely 7K in T^ff and 0.15% in bolometric flux have been 
found, consistent with the effect that the absolute calibration in 
iRI)c can introduce. For the sake of homogenizing the stellar pa- 
rameters derived in this work, also these small offsets have been 
corrected for stars with no iRI)c photometry. 

The effective temperature can be determined from any in- 
frared photometric band, in our case JHKs from 2MASS. 
Ideally all bands should return the same T^ff, but photometric 
errors and zero point uncertainties in the absolute calibration of 
each band introduce random plus systematic differences. In the 
case of 2MASS, those amount to few tens of K as we show later 

The magnitude in a given band ^ is converted into a physical 
flux (i.e. erg cm"^ s"' A"') via 



!rf(Earth) = 57"'(Earth)10-"'*<"f-<' 



(1) 



which depends on the zero point (m^"*) and the absolute flux cal- 
ibration {7^^''^) of the standard star defining the photometric sys- 
tem under us^ 

Most of the photometric systems, including Johnson- 
Cousins and 2MASS, use Vega as the zero point standard. Vega's 
flux and magnitudes in different bands have been notoriously 
difficult to measure with sufficient accuracy (e.g. iGray 2003 
and references therein). The problem is only apparently resolved 
when resorting to Robs'- in the ideal case of a unique template 
spectrum for Vega the choice of its absolute calibration would 
cancel out in the ratio. In practice, the situation is far from this 
since the pole-on and rapidly rotating nature of this star imposes 
the use of a composite abso lute cahbrated spe ctrum for differ- 
ent wavelength regions (e.g. ICasagrande et al.ll2006. and refer- 
ences therein). Such complication does not disqualify Vega as a 
spectrophotometric standard, but it makes its use more problem- 
atic. From Eq. ([T) it can be immediately noticed that a change 
of 0.01 mag. corresponds to a change of about 1% in flux. Since 
it is possible to interchangeably operate on both zero points and 
fluxes, for the sake of our discussion it is their composite effect 
that must be considered, though in the following we shall usually 
refer to fluxes. 

Recently, HST spectrophotometry for Vega has provided a 
unique calibrated spectrum extending from 3200 to 10000 A 
with 1 - 2% accuracy (Bohlin 2007). In the infrare d, once the 
zero points newly determined from Maiz-Apellaniz (2007) are 
used, this result is also in broa d agreement with the 2MASS ab- 
solute calibration provided by ICohen et al] (|2003|) . iRieke et al] 
(l2008h have also recently reviewed the absolute physical cali- 
bration in the infrared, substantially validating the accuracy of 
2MASS: their recommended 2% increase of flux in Ks band is 
in fact compensated by their newly determined zero point for 
Vega, thus implying an effective change in the overall Ks cali- 
bration of only 0.2%. We have tested all these different possibil- 
ities; wi th respect to the HST and 2MASS calibration adopted in 
ICasagrande et al.i (l2006i) the derived T^ff are affected at most by 
20 K. Such difference is thus within the aforementioned global 
2% uncertainty which allows for systematics in Tf^g of order 



' We point out that Eq. l[T) holds exactly for a heterochromatic mea- 
surement, while for computing a monochromatic flux from the observed 
photometry, an additional correction (the so called ^-factor) must be 
introduced to account for the fact that the zero point of the photomet- 
ric system is defined by a standard star, which usually has a different 
spectral energy di stribution across the filter window with resp ect to the 
problem star (e.g. lAlonso et al.l[l996al : ICasagrande et alJ2006h . 
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40 K. Our zero points an d absolute fluxes a re ess entially iden- 
tical to those adopted in ICasagrande et alJ (l2006l) except for a 
small fine-tuning which will be further discussed in Section[3] 

Despite the recent increasing concordance in establishing ab- 
solute fluxes, the uncertainties which have historically plagued 
Vega are crucial in the context of understanding the effective 
temperatures determined via IRFM by various authors. We have 
tested that uncorrected changes of a few percent in the absolute 
calibration of optical bands (needed to recover the bolometric 
flux) can introduce spurious trends with Teff and [Fe/H] up to 
few tens of K. Similar changes in the absolute calibration of in- 
frared bands have only minor impact on the bolometric flux, but 
as akeady mentioned, Teff is very sensitive to them since they 
enter explicitly in the definition of Robs' increasing all of them 
by 2% translates into a decrease of approximately 40 K in Teff. 
Considering that differences of few percent in the adopted zero 
points and fluxes are commonly present among various IRFM 
implementations, it can be immediately realized that they are re- 
sponsible for systematic differences among various authors. 



2.3. Alonso etal. (1996) scale 

One of the most ext ensive applica t ions of the IRFM to Pop I and 
II dwarfs is that of I Alonso et al ] (Il996ah . which was based on 
the infrared photometry collected at the TCS (Telescopio Carlos 
Sanchez, Alonso et al. 1994b) and absolutely calibrated using 
a semi-empirical approach relying on (mostly) giant stars with 
measur ed angular diameters to determine the reference absolute 
fluxes dAlonso et al.l 1994ah. The comparison between our T^efF 
and those bv lAlonso etal Jo 996a') is shown in Figure[3] Despite 
the scatter arising from the different input data we used, there 
is a clear offset with our scale being systematically hotter. No 
obvious trends in T^ff and [Fe/H] appear. This offset is easily 
explained in terms of the absolute calibration underlying the two 
different photometric systems adopted. This involves the trans- 
fo rmation from TCS to 2M ASS system (see also the discussion 
in ICasagrande et a"l]|2006h . which could in principle introduce 
additional noise (see Section 12741 ). A more detailed description 
of the absolute calibration (and angular diameters) employed by 
Alonso and a comparison with our own is presented in Appendix 
A. 



An area of particular interest is the determination of effective 
temperatures in very metal -poor, turn-off star s . We h ave tested 
the effect of using the new Castelli & Kuruczl ( 2004b model at- 
mospheres in the I RFM instead of the lKunicz ([1993) adopted by 
lAlonso et all ('1996a). The IRFM is known to be little model de- 
pendent (e.g. Asplund & Garcia Perez 2001; Casagrande 200^ 
and in fact there a re no big differences exc ept at the lowest 
metallicities, where ICastelli & Kuruc j (|2004|) support effective 
temperatures hotter by ~ 40 K. The reason for such a dis- 
crepancy stems from the new models returning higher flux be- 
low ~ 4000 A, a region where the most metal-poor, turn-off 
stars commence emitting non negligible amounts of energy. 
Since we do not have UV photometry (and its standardiza- 
tion would be uncertain), we must rely on model atmospheres 
to determine the flux over this region (Figure HJi. The latest 
model atmosphere calculations show excellent agreement as we 
checked that near ly identical T^g are ob tained when the new 
MARCS models (Gustafsson et al. 2008) are used instead of 
those b y Castefli & Kurucz (2004) (also Section |53?T] ). but see 
lEdvardss on (2008) for a discussion of the performance of model 
atmospheres in the blue and ultraviolet. 
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Fig. 3. Difference between the effectiv e temperatures obtained in 
this work (TW) and those reported in lAlonso et al. I ([T996ah for 
220 stars in common. In case of reddening, only stars with val- 
ues of E{B - V) equal to within 0.02 mag. have been plotted. 
Thick continuous lines connect the means computed in equally 
spaced bins of [Fe/H] and Teff. Error b ars are t h e stan dard de- 
viation in each bin. Top panels: when 'Kuruc^ (Il993h models 
are used in our version of th e IRFM. Bottom panels: when the 
new ICastelli & Kuruc j (2004) models are used instead. Below 
[Fe/H] = -1.5 the new models support T^ff hotter by 20 to 40 K. 



2.4. Ramirez & Melendez (2005) scale 



A revision of the lAlonso etal. (Il996a|) implem e ntation of the 
IRFM was carried out by Ramirez & Melendej (l2005a) based 
on the TCS (for the computation of Rtheo) and Johnson's (for 
the c omputation of the bolometric fluxes) JHK photometric sys- 
tems (lAlonso et al.ll994 b':'Bessen & Brett' 1988). Here we repli- 
cate the TefF determination by Ramirez & Melendez (2005a) for 
comparison purposes. When running their implementation, we 
transformed the 2MASS photometry into TCS using their equa- 
tions. However, when comparing the transformed and original 
JHK values for these stars we found zero point differences at the 
level of 0.01 magnitudes: these offsets are within the photometric 
uncertainties and smaller than the scatter in the fits leading to the 
transformation equations, but they introduce changes in the de- 
rived Teff values up to few tens of K (see Section l272l i. Therefore 
we took those into account to precisely transform 2MASS data 

into the TCS system. 

The iRamfrez & MelendezI (12005 ah bolometric fluxes were 
determined using the K-band bolometric correction calibration 
bv lAlons o et al. (1995), which depends only on the Johnson 
(V - K) colour index and the stellar metallicityQ This cali- 
bration is internally accurate within its ranges of applicability 
and one would expect that extrapolations slightly outside these 
ranges would still provide reliable results at low metallicities. 
This approach was followed by Rarm'rez & Melendez (2005^. 
With regards to the a bsolute flux calibration in the infrared, 
Ramirez & Melende"3 (l2005al) adopted that of lAlonso et"an 
( 1994a), which is vaUd for TCS JHK photometry while we use 



We have also tested that in the context of computing 
bolometric fluxes for this work, the updated J. Carpenter 
transformations from 2MASS to Johnson available online at: 
http://www.astro.caltech.edu/~jmc/2mass/v3/transformations are 
instead accurate enough and insensitive to small zero point changes. 
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Teff=6000 K, [Fe/H] = -4.0 Teff=6000 K, [Fe/H] = -2.5 Teff=6000 K. [Fe/H] = 0.0 
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Fig. 4. Comparison between fKurucj (Il993h (thick line) and'Castelli & Kurucd (|2004') (thin line) synthetic spectra at different metal- 
licities for an assumed log g = 4.0. Shaded area is the wavelength region covered by our multiband photometry. The difference in 
the UV flux gets more prominent when going to more metal-poor stars, but for the sake of the IRFM is entirely negligible at solar 
metallicity. 
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Fig. 5. Top panels: difference between the effective temperatures 
of this work (TW) and those obtained when the same input 
data are used in the Ramirez & Melendez (2005a) implemen- 
tation (RM05). Bottom panel s: as in the top panels but for the 
iRamirez & Melendej (12005 al ) temperatures re-determined using 
the bolometric fluxes obtained in this work. 



an update of ICohen et aP (l2003h for the JHKs 2MASS system 

(see also Section|3]l. 

The difference between our results an dH amfrez & Melende"3 
(l2005ah when the same input data and reddening values are 
adopted is illustrated in the top panels of Fig.|5] Some of the scat- 
ter arise from transforming 2MASS magnitudes into TCS, but 
clear trends with both with Tetf and [Fe/H] are present. For the 
bulk of the stars with [Fe/H] > -2.0 and 4800 < T^ff < 6200 K a 
roughly constant offset of about lOOK is observed, our stars be- 
ing hotter. In the metal-rich regime such an offset is present also 
for hotter stars (Teff > 6200 K), but reduces somewhat for the 
coolest metal-rich dwarfs, reaching a minimum of about 50 K 
at Teff - 4500 K. A steep trend is seen for moderately metal- 
poor dwarfs (-2.0 < [Fe/H] < - 1.0) below 4800 K a region 
with few or no calibrating stars in lAlonso et al ] (ll995h . For the 
warmer, most metal-poor stars in the sample, the differences 
decrease sharply with increasing Teff and decreasing [Fe/H], 
quickly becoming negative i.e., Ramirez & Melendez (2005a) 
temperatures become warmer, reaching a maximum value of 
about -lOOK at T^ff ^ 6500 K and [Fe/H] -3.5. 



To investigate the source of these differences, we re- 
calculated the IRFM temperatures of Ramirez & Melende^ 
(2005a) using our b olometric fluxes instead of th e calibration 
formulae adopted bv lRamfrez &Melendezl(l2005ah . This choice 
is perfectly legitimate, since what is crucial in the IRFM are 
the infrared fluxes which appear explicitly in the definition of 
Robs, while Teff depends only mildly on the bolometric flux 
(Section [2. 2l i. Therefore, adopting our bolometric fluxes is sub- 
stan tially independent of the und erlying temperature scale, i.e. 
the Ramirez & Melendezl(l2005ah scale is still recovered despite 
now using the new bolometric fluxes determined in the present 
work. The result of this exercise is shown in the bottom pan- 
els of Fig. |5l The m ajor trends caused from extrapolating the 
lAlonso et all (fT99l bolometric formulae now disappear with a 
constant offset AT^s = 85 + 13 K above 5000 K. The smafl trend 
that rema ins below this temperature c orresponds to the thresh- 
old where lRamfrez & Melende3(l2005 a) stop using the / band to 
determine T^s, which in the TCS system usually returns slightly 
cooler Teif than H and K bands. 

Fro m this comparison it is clear that iRamfrez & Melendez! 
(l2005ah temperatures for the metal-poor turn-off stars are 
warmer due to the use of a photometric calibrat ion to derive 
the bolometric fluxes. In fact, we realize that the lAlonso et al.l 
(Il99 5) formula is robust down to [Fe/H] ^ -2.5 and up to 
Tgff - 6500 K but only a few c alibrating stars more metal-po or or 
warmer exist in their sample. iRamirez & Melendez! (l2005al) use 
of this formula in regions where the calibration is uncertain (and 
in some cases outside of the ranges of applicability) has resulted 
in the very high temperatures of the more metal-poor turn-off 
stars. The extrapolation is, of course, not a valid procedure, even 
though one might expect the [Fe/H] dependence of the calibra- 
tion not to be so important at low metallicity. However, as can be 
seen from Fig. 4 in Alonso et al. ( 1995), at these relatively high 
temperatures, the effect of [Fe/H] is very important and such 
extrapolations should not be performed. 

The difference that remains after adopting consistent bolo- 
metric fluxes between this work and Ramirez & MelendeJ 
(12005 a!) (lower panels of Fig. |5]l is mostly due to the use of dif- 
ferent infrared absolute flux calibrations. In fact, by lower ing 
the absolute fluxes adopted bv lRamfrez &Melendezl ( l2005ah bv 
about 4%, the mean difference reduces to almost zero. We thus 
conclude that our and Ramirez & Melendez (2005a) IRFM im- 
plementations can be made perfectly compatible if the same in- 
put parameters and flux calibration are used. 
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Fig. 6. Top panels: difference between the effec- 
tive temperatures of this w ork (TW) and those in 
iGonzalez Hernandez & Bonifaciol (|B009) (GB09) for 380 
stars in common. Filled circles are stars with Teff < 5000 K 
without {RI)c photometry in GB09. Bottom panels: as in the top 
panels, but when the same reddening corrections are used. 



2.5. Gonzalez Hernandez & Bonifacio (2009) scale 

The most recent work on the IRFM is that by 
iGonzalez Hernandez & Bonifaciol (l2009l) . which is also based 
on 2MASS photometry. The main difference between theirs 
and our implementation is the different absolute calibration and 
zero points adopted for Vega. They based their work on the 
[Castelli & Kurucz ( 1994) model and McCall (2004) magr iitudes 
instead of th e HST (LB ohhn & Gilliland 120041: iBohlinI |2007|) 
and 2MASS (ICohen et al...2003.) values that we use. Although 
such differences are within the current observational errors, in 
the infrared the combined effect of their fluxes and zero points 
is on average 1.5 - 2.0% higher than ours, implying effective 
temperatures cooler by 30 - 40 K (see Appendix A). This can 
be immediately appreciated in Figure |6] which indeed shows a 
constant offset of this magnit ude for stars in common, thus con- 
firm ing the offset noticed by Gonzalez Hernandez & Bo nifaciol 
(I2009.) for stars in common with Casagrande et al. (2006i) . 

The very steep trend at the lowest metallicities is due to the 
different reddening corrections we adopt with respect to theirs. 
When the same E(B - V) values are adopted (bottom panels 
in Figure |6]l, the offset remains constant throughout the entire 
[Fe/H] and Tetf range, except for few outliers due to the differ- 
ent input data (mostly optical photometry) adopted. This clearly 
stresses the importance of proper reddening correction for de- 
termining effective temperatures via IRFM in stars outside of 
the local bubble. For the most metal-poor stars in the sam- 
ple, we use interstellar NaD lines to achieve higher precision 
(Section lZTb while Gonzalez Hernandez & Bonifacio (2009) re- 
sorted to reddening maps scaled by the distance and the galac- 
tic latitude of the star and scale height of the dust layer. The 
trend towards cooler effective temperatures that we obtain in 
this regime thus stem entirely from better reddening corrections. 
Finally, we suspect that the tr end for T^tf < 5000 K is due to the 
absenc e of iRI)c colours in IGonzalez Hernandez & Bonifaciol 
(l2009l) (SectionO bottom panel of Figure|2ll. 



3. Resolving different versions 

It is clear from the discussion above that we now understand 
where different T^ff scales originate from and the crucial role 
played by the absolute calibration. Our approach has been to 
adopt the latest calibration available for each photometric sys- 
tem: currently those are accurate at the 2% level, implying pos- 
sible systematic uncertainties of order 40 K. Here we want to 
improve upon this uncertainty using an independent verification 
of the absolute calibration adopted. 



3.1. Solar twins 

The use of solar-type stars to calibrate photometric systems 
has a long and noble history, which relies on taking absolutely 
calibrated measurements of the Sun and computing s ynthetic 
colours to compare with other solar -type stars (e.g. iJohnsonI 
ll965l:ICampins et all 1 9851: iRieke et al.i.2008) . This rationale can 
be extended to other physical properties, namely using the solar 
effective temperatur e Tss.c = 5777 K a s the average value for 
solar-type stars (e.g. lMasana et al.ll2006l) . This technique is well 
established and goes under the name of solar analogs method, 
but there is some sort of petitio principii in the underlying Teft 
scale adopted and/or the solar colours assumed to select solar 
analogs in first instance. 

A way to break such a degeneracy is provided by solar 
twins, i.e. stars with spectra indistinguishable from the Sun 



dCavrel de Strobel & Bentomalll989l:|Porto de Mello & da Silval 



119971) . Our twins were drawn from an initial sample of about 
100 stars broadly selected to be solar like: the identification of 
the best ones was based on a strictly differential analysis of high- 
resolution (R ~ 60000) and high signal-to-noise (S /N > 150) 
spectra with respect to the solar one reflected from an aster- 
oid and observed with the same instrument. Within this initial 
sample, the selection criterion adopted to identify the best twins 
did not assume any a priori effective temperature or colour, 
but was based on the measured relative difference in equiva- 
lent widths and equivalent widths vs. excitation potential re- 
lations with respect to the observed s olar reference spectrurn 
and thus entirely model in dependent dMelendez et alj l2006al: 
iMelendez & Ramfrezll2007l) . Since the spectra of the solar twins 
match so closely the solar one, exceedingly accurate differential 
spectroscopic analysis with respect to Teff o, [Fe/H]^ and log^o 
is possible (Melendez et al. 2009a; R amirez et al . 200^. 

Ten stars were identified as most closely resembling the Sun 
and are given in Table|2] including HIP56948, the best solar twin 
curre ntly known (Melendez & Ramirez 2007; Takeda & Tajits^ 
l2009b . A crucial requirement for these stars is to have accu- 
rate and homogeneous photometry in order to derive reliable 
Teff via IRFM. While this is possible in the infrared because of 
2MASSQ, optical photometry is also important to properly re- 
cover the bolometric flux where these stars emit most of their 
energy. Johnson-Cousins photometry would be the ideal choice, 
but unfortunately is not available for all these targets. To over- 
come this limitation, in the optical we used the Tycho2 BtVt 
system which uniformly and precisely covers the entire sky in 
the magnitude range of our interest (H0g et al. 2000). Notice that 
we did not transform BtVt into BV but instead implemented 
our IRFM code to work directly on the Tycho2 system. Also, 
as discussed in Section \22] the absence of iRI)c photometry is 
not relevant for stars hotter than 5000 K. All twins are closer 



•* In fact, the other we ll known solar twin 18 Sco 
iPorto de Mello & da Silvalll997D has saturated 2MASS colours. 
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than 72 pc, where reddening is expected to be zero or negligi- 
ble: nearly all of them have Stromgren photometry (Melendez 
et al. in prep.) and the Schuster & Nissen ( 198§J reddening cah- 
bration confirms indeed such a conclusion. 

3.2. A finely tuned absolute calibration 

As for the Johnson-Cousins system, we based the absolute cali- 
bration ofdie Tycho2 system on Vega (iBohlin & Gilliland 2004; 
[Bohlin 2007"), adopting the BtVt zero points of iMafz- ApellanizI 

(2007) and the corresponding filter transmission curves of 
lBesselll(l2000l) . 

The first instance, we determined Teif via IRFM for each of 
the twins in Table |2l their average eff'ective temperature turned 
out to be 5782 K, remarkably close to Tetf o, thus confirming the 
high accuracy achieved using the HST and 2MASS absolute cal- 
ibration. Based on Monte Carlo simulations with the photomet- 
ric eiTors in Table|2] the uncertainty in Teff determined via IRFM 
is of order 30 K for single stars. Imposing the mean effective 
temperature of all solar twins to equal Teff o we estimate the un- 
certainty on the zero point of our temperature scale to be 15 K 
based on a bootstrap procedure with one million re-samples. At 
the same time, for HIP56948 we also recover Teff ^ within Icr. 

Though the solar twins test confirms the global reliability 
of the adopted absolute calibration, for all stars in Section ITTI 
having Tycho2 photometry and Teif > 5000 K we further re- 
quired each infrared band to return on average the same Tetf as 
the others (Figure |7]i. By imposing such a consistency we im- 
prove upon small systematic trends which could arise when de- 
termining effective temperatures in stars with Teff and [Fe/H] 
very different from our solar twins. This led to a decrease of 
the absolute calibration by 1.6% in the J band and an increase 
by 1.5 and 0.3% in the H and Ks bands, respectively (see also 
Appendix A). In terms of synthetic magnitudes these differences 
make H and Ks redder by 0.016 and 0.003 and J bluer by 0.017, 
th us removing alm ost entirely the infrared colour offsets found 
bv lCasagrande et a l. (2006) when comparing observed and syn- 
thetic photometry. We cannot entirely rule out whether these 
systematic differences arise from the adopted synthetic library 
or the absolute calibration, but since the IRFM depends only 
marginally on model atmospheres and the infrared spectral re- 
gion is relatively easy to model, we are strongly in favour of 
the second possibility. From a pragmatic point of view, this im- 
proves the consistency in determining Tgff. Also, such changes 
are within the 2MASS quoted errors and for the Ks band we re- 
mark the agreement with the 0.2% increase found bv lRieke et al.l 

(2008) and discussed in Section IZ2] As expected, stars with the 
best 2MASS pedigree also return better agreement in all bands 
(full circles in Figure|7]l. We have also checked that the increas- 
ing scatter in Figure|7]is prim arily due to photometric eiTors. We 
recall that iRieke et all 12008 ) found a 2% offset between Read 
1 and Read 2 mode in 2MASS|fl though they were not able to 
derive a universal correction for this effect. All our solar twins 
have Read 1 mode and the absence of a universal correction sug- 
gests that while Read mode 2 can decrease the precision of Teff 
the overall accuracy of our calibration remains valid. 

With the fine-tuning discussed above, the median (mean) 
effective temperature of our solar twins is 5777 (5779) K. 
Restricting only to the twins having Teff within the observa- 
tional errors, still confirm such conclusion. As a further indepen- 
dent test, we applied our IRFM to the list of solar analogs used 



^ This mode indicates which readout is used to derive photometry 
|http://www.ipac.calt ech.edu/2mass/releases/allsky/doc/sec3-lb. html, 



bv lRieke et al.l (l2008h and determined their median (mean) T^ff 
to be 5791 (5786) K, thus confirming the reliability of the zero 
point of our temperature scale, which has an uncertainty of 15 K. 
Such a value implies possible systematics in the absolute cali- 
bration at the 1 % level. The systematic error in recovering the 
bolometric luminosity is however smaller since infrared fluxes 
enter twice in Rohs, thus partly compensating their uncertainty. 

The corrections in the infrared absolute calibration discussed 
here have been used also in determining T^ff for stars in Section 
12.11 Since for those stars we are using Johnson-Cousins pho- 
tometry, there could still be small differences arising from the 
absolute calibration in the optical: for stars in common a mean 
systematic of 8 K in T^s and 0.15% in bolometric flux was found 
and corrected. 

4. Validating the proposed temperature scale 

The IRFM determines T^ff in an almost model independent way, 
primarily recovering the bolometric flux ;FBo/(Earth) of the star 
under investigation. From the basic definition linking those two 
quantities the stellar angular diameter 0irfm can be obtained self- 
consistently and this was a ctually one of the driving re asons for 
developing the technique (i Blackwell & ShailisllI977 '). In what 
follows, we use this information to further validate our results. 

4. 1 . Interferometric angular diameters 

An independent test of accuracy for the zero point of our effec- 
tive temperature scale involves the comparison with the angular 
diameters measured using interferometric techniques (coiTected 
for limb-darkening, hereafter denoted by 0ld)- In our case, angu- 
lar diameters are a natural consequence of the Teff determination 
procedure and for each star the Teff , T'boI, ^irfm values are self- 
consistent, i.e., they represent a unique solution for a given set of 
input data. We also prefer to compare angular diameters directly 
(i.e. ^iRFM vs. 0ld) since the effective temperatures reported in 
various interferometric works would be more heterogeneous be- 
cause of the adopted bolometric coiTections. 

Given the difficulties involved in the measurement of the 
small angular diameters of dwarfs and subgiants (even the near- 
est ones have angular diameters below 10 milli-arcseconds), 
only a relatively small group of such stars has been observed to 
date for that purpose (see also Appendix A for a discussion of the 
angular diameters used by Alonso et al. 1994a). We performed 
a literature search for interferometrically determined angular di- 
ameters with precision better than 5 % (which corresponds to an 
accuracy of 2.5% in effective temperatures, roughly 150 K at so- 
lar temperature, assuming no error in the bolometric flux) and 
found data for 28 stars, 16 of which have 0ld measured to bet- 
ter than 2 % (Table [3]i. The efforts made by the interferometry 
community in the last few years are commendable given that the 
nu mber of stars with reliable Old has nearly doubled since 2005 
(cf. lRamfrez & Melendejl2005ah . 

Unfortunately, all dwarfs and subgiants with reliable Old are 
brighter than V - 6, implying infrared magnitudes < 5 where 
2MASSphotometry has large observational eiTors and starts to 
saturatqj Therefore we cannot apply our IRFM directly on them 
to get 0IRFM- Instead, we adopt an indirect approach using the 
photometric reff:colour and fsohcolom relations presented in 
Section |6] Using the photometry of our sample stars (i.e. those 
used in the construction of the calibrations and therefore with 
Teif directly determined via IRFM), we checked that the zero 

* www.ipac.caltech.edu/2mass/releases/allsky/doc/sec2_2.html#pscphotprop 
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Table 2. Tycho2 and 2MASS photometry for our solar twins sample, together with the spectroscopic parameters and the effective 
temperatures determined via IRFM. For the latter, the errors are those arising from the photometry alone, not including the 15 K 
uncertainty in the zero point of our temperature scale. All twins have "A" quality flag and Read 1 mode in all 2MASS bands. 
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Fig. 7. Top panels: difference between T^ff and the effective temperature determined in each infrared band before tuning the absolute 
calibration. Full circles are stars with quality flag "A", Read 1 mode and total 2MASS photometric errors < 0.07 mag. while open 
circles are for all other stars. Bottom panels: as in the top panels, but with the adjusted absolute calibration. H band photometry has 
usually slightly higher error than J and Ks and the final temperature is the weighted average of that obtained in each band. 



point of our Teff and TboI scales is correctly reproduced by the 
calibration formulae presented in Section|6] independently of the 
apparent magnitudes of the stars. Also, for the two stars having 
HST spectrophotometry (next Section) we checked that our cal- 
ibration formulae reproduce nearly the same results as directly 
applying the IRFM. We were careful about propagating all pos- 
sible sources of random error such as uncertainties in the input 
photometry, metallicity, and the reliability of the colour calibra- 
tions, as quantified by the standard deviation of each polyno- 
mial fit (Tables|4]and|5]l. For most of the stars with reliable 0ld 
(i.e. better than 2%), only BV photometry was available, while 
for the remaining BViRFfc was u sed. Metallicities were adopted 
from the updated version of the ICavrel de Strobel et al ] (1200 Ih 
[Fe/H] catalog by Melendez (in prep.), which nearly triples the 
number of entries in the original catalog. 

The comparison of the angular diameters measured inter- 
ferometrically with those derived using our IRFM colour cali- 
brations is shown in Fig. |8] (see also Table O. Stars that have 
0LD determined with accuracy better than 2 % are shown with 
full symbols. Using only the latter, the average difference in 
angular diameter (IRFM-LD) is -0.62 + 1.70% which corre- 
sponds to a zero point difference in the effective temperature 
scale of only +18 + 50K at solar temperature. This is also in 
agreement with the uncertainty on the zero point of our tem- 
perature scale discussed in Section 13.21 No obvious trends are 



seen with [Fe/H] (from about -0.8 to -hO.3) or T^s (from 4400 
to 6600 K). Note, however, that if we exclude the two coolest 
stars (from the group of those having errors smaller than 2 %), 
a small trend is seen with Teff. The trend -if real- appears more 
clearly for early type stars, with ^irfm being underestimated (and 
therefore the IRFM effective temperatures overestimated) with 
respect to the interferometric measurements. Interferometry re- 
sorts on ID model atmospheres to correct from the measured 
uniform-disk angular diameter to the physical limb-darkened 
disk to which we compare with. Interestingly, 3D models pre- 
dict less center-to-limb variation tha n ID m odels as moving 
from K to F type stars (Allende Priet o et al] 12002; Big ot et al] 
'2006). Reduced limb-darkening corrections imply smaller 6u): 
the trend discussed above qualitatively fit into this picture. How 
well our result agrees quantitatively with this picture we leave to 
future studies. 

Interestingly, 'Rami rez & Melendez! (l2005al) made a similar 
comparison of angular diameters and also found good agree- 
ment with their IRFM T^s scale, which is, however, systemat- 
ically co oler (by - 1 00 K) than the present one for [Fe /H] > -2 
(see also ICasagrandel2008l) . We compared the stars with angular 
diameters in common between table 4 of Ramirez & Melendez! 
(2005a, RM05) and the present study (C09, TableO and found 
an average difference (C09-RM05) of 0.1 + 2.2% in angular 
diameters, 3.0 + 3.0% in bolometric fluxes and 40 + 37 K 
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iRamfrez & Melended (l2005ah (see the discussion on the small 
zero point differences to convert 2MASS into TCS presented 
in Section IZ4l i. To g auge further insights, w e redetermined the 
temperatures used by iRamfrez & MefindezI (l2005al) using their 
colour calibrations for the same BV{RI)c input data we adopted 
in this section and found AT^ff - 72 + 52 K. In addition, 
we adopted our bolometric fluxes lowered by 1 %, which cor- 
responds to the average difference we find for our complete 
sample. In this case the difference in angular diameters sets to 
-2.4 + 2.1 %, much closer to the expected -3 %, offsetting the 
iRanurez & Me lendez (2005a) scale with respect to interferomet- 
ric measur ements. S ince the present work represents an improve- 
ment over IRamfrez & Melen dez (2005a), in particular the fact 
that the T^ff, 'FboI, ^rfm values are a self-consistent and unique 
solution to each problem star, and given that the number of com- 
parison stars has doubled since 2005 (note also that the 0ld val- 
ues of some stars have been re-determined), it is likely that the 
good agreement found bv .Ramfrez & Melendez (.2005a) was due 
to a conspiracy of photometric errors which propagated to both 
Teff and TboI determinations and low number statistics. More 
measurements of stellar angular diameters via interferometry 
are clearly necessary, and therefore highly encouraged, to bet- 
ter constrain indirectly determined effective temperature scales. 
However, as this exercise has shown, many critical ingredients 
enter in the comparison with angular diameters. In particular 
bolometric corrections and effective temperatures should be de- 
termined as self-consistently as possible, also avoiding transfor- 
mation between photometric systems. It gives us confidence that 
the zero point uncertainty from solar twins, angular diameters 
and HST spectrophotometry (next Section) returns in all cases 
independent and very consistent results. 

While the angular diameter comparison does not extend be- 
low [Fe/H] ^ -1.0, leaving our results for halo stars "un-tested" 
in this context, in the next Section we use HST spectrophotome- 
try to gauge further insight on the topic. 
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Fig. 8. Top two panels: Comparison of angular diameters mea- 
sured interferometrically (0ld) and via our IRFM photometric 
calibrations (0irfm)- Full symbols represent stars that have 0ld 
measured with accuracy better than 2%. Bottom two panels: 
Difference (in %) between 0ld and 0irfm as a function of stellar 
parameters. Solid lines represent 1-to-l correspondence, dashed 
and dotted lines are the average difference and 1-cr error for the 
full data point, respectively. 



in Teff. Given the large scatter, these numbers are still con- 
sistent with the mean differences in T^s and TboI from these 
two studies (Section 12.4b . however, we would expect our di- 
ameters to be roughly smaller by 3 %, our fluxes brighter by 
1 % a nd our Teff hotter by 100 K (see also ICasagrande et alJ 
l2006l) . While TboI and T^s compensate to give almost exactly 
the same angular diameters, the 40 K offset might be more rep- 
resentative of the difference with the TCS magnitudes used in 



4.2. HST spectrophotometry 

For each star, we obtain a synthetic spectrum tailored at the ef- 
fective temperature determined via IRFM (Section 12.21 ). Since 
the angular diameter is determined, each synthetic spectrum is 
absolutely calibrated (i.e. in units of ergcm^^s"' A"'), and can 
be used to further test our results. In fact, from F- to early K-type 
stars, all continuum characteristics approximately longward of 
the Paschen discontinuity depend almost exclusively on the ef- 
fective temperature, relatively unaffected by spectral lines and 
NLTE effects as well as from the treatment of convection. 

The CALSPECQ library contains composite stellar spectra 
measured by the STIS (0.3-1 .0 fim) and NICMOS ( 1 .0-2.5 yum) 
instruments on board of the HST and used as fundamental flux 
standard. Free of any atmospheric contamination the HST thus 
provides the best possible spectrophotometry to date, with 1-2% 
accuracy, extending from the far-UV to the near infrared. The 
absolute flux calibration is tied to the three hot, pure hydrogen 
white dwarfs, which constitute the HST primary calibrators, nor- 
maHzed to the absolute flux of Vega at 5556 A(Bohlin 2007|). 
Thus, except for the normalization at 5556 A the absolute fluxes 
measured by STIS and NICMOS are entirely independent on 
possible issues regarding Vega's absolute calibration in the in- 
frared and offer an alternative approach to the 2MASS calibra- 
tion underlying our temperature scale. 



' http://www.stsci.edu/hst/observatory/cdbs/calspec.html 
January 2009. 
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Fig. 9. Left upper panel: comparison between the observed HD209458 CALSPEC spectrum (black line) and the synthetic spectra 
derived for two different Teff , using our preferred absolute calibration (blue line) and increasing the infrared absolute calibration by 
5% (red line). Left lower panel: ratio of synthetic to observed spectra. Full circles are the ratio between the fluxes obtained once 
the Vega calibration is used with the observed magnitudes and the fluxes obtained directly from the convolution of the CALSPEC 
spectrum with the appropriate filter transmission curve. Error bars take into account uncertainty in the Vega calibration and zero 
points, as well as in the observed magnitudes. Right panel: reduced for various Tetf solutions corresponding to different adopted 
absolute calibrations. Our choice (Section |372] | always lies very close to the minima obtained fitting a parabola to the data (lines of 
different style). Different symbols correspond to cut longward of 0.66 //m (diamonds), 0.82 /im (squares) and 1.46//m (triangles) 
as explained in the text. The sigma levels have been computed using the incomplete gamma function for the number of degrees of 
freedom longward of our cuts. 



Table 3. Stars with measured interferometric angular diameters. 



are available: th e exoplanet host star HD209458 (e.g 

ICharbonneau et all 12000) and the fundamental SPSS stan dard 



HD 




Ref." 


cff 


[Fe/H] 






mas 




K 


dex 


mas 


3651 


0.790 ± 0.027 


1 


5234 


0.15 


0.756 ± 0.022 


6582 


0.973 + 0.009 


2 


5403 


-0.84 


0.954 ±0.021 


9826 


1.114 + 0.009 


1 


6151 


0.10 


1.121 ±0.023 


10700 


2.078 ±0.031 


3 


5364 


-0.53 


2.089 ± 0.026 


10780 


0.763 ±0.021 


2 


5317 


0.01 


0.806 ± 0.022 


19994 


0.788 ± 0.026 


1 


6020 


0.18 


0.746 ± 0.009 


22049 


2.148 ±0.029 


3 


5056 


-0.09 


2.200 ± 0.032 


23249 


2.394 ± 0.029 


3 


5060 


0.08 


2.399 ± 0.059 


26965 


1.650 ±0.060 


3 


5188 


-0.27 


1.482 ±0.018 


61421 


5.443 ±0.030 


3 


6626 


0.00 


5.326 ±0.068 


75732 


0.854 ± 0.024 


1 


5282 


0.38 


0.718 ±0.025 


102870 


1.450 ±0.018 


4 


6100 


0.13 


1.426 ±0.014 


117176 


1.009 ±0.024 


1 


5540 


-0.06 


0.969 ±0.021 


120136 


0.786 ±0.016 


1 


6407 


0.28 


0.840 ± 0.019 


121370 


2.244 ±0.019 


3 


6052 


0.26 


2.214 ±0.043 


128620 


8.511 ±0.020 


3 


5772 


0.20 


8.511 ±0.079 


128621 


6.000 ±0.021 


5 


5217 


0.23 


6.151 ±0.234 


131977 


1.230 ±0.030 


3 


4633 


0.04 


1.162 ±0.054 


150680 


2.397 ± 0.044 


3 


5780 


0.03 


2.352 ± 0.055 


161797 


1.953 ±0.039 


3 


5520 


0.22 


2.004 ± 0.050 


185144 


1.254 ±0.012 


2 


5293 


-0.21 


1.261 ±0.029 


188512 


2.180 ±0.090 


6 


5164 


-0.18 


2.070 ± 0.049 


190360 


0.698 ±0.019 


1 


5564 


0.21 


0.673 ± 0.017 


198149 


2.650 ± 0.040 


3 


4980 


-0.16 


2.720 ± 0.090 


201091 


1.775 ±0.013 


3 


4429 


-0.24 


1.706 ±0.070 


209100 


1.890 ±0.020 


3 


4665 


-0.06 


1.825 ±0.021 


217014 


0.748 ± 0.027 


1 


5754 


0.17 


0.698 ±0.019 



1.- [ Baines et al ^ ('2008'), 2.- iBovaiian et all j2008h . 3.- 
iKervell a & Fouque ( 200^ (weighted average if more than 
one measuremen t wa s available), 4 . - iNorth et alJ ( l2009h . 5.- 
iBigot et al.l j2006h . 6.- lNordgren et alj ( 119991) . 



Two of the CALSPEC targets are late-type main-sequence 
dwarfs for which accurate photometry, \ogg and [Fe/H] 



BD H-174708 (e.g. Fukugitaet al. 1996; Smith etaDllQOl. For 
each of these targets we computed T^a and derived the corre- 
sponding physical flux using the absolute calibration presented 
in Section [3T2l For comparison, we also determined the effective 
temperatures and the corresponding fluxes when changing our 
adopted infrared absolute calibration by different amounts up to 
+5%, which roughly correspond to +100K in T^g. The agree- 
ment was quantified using statistics between the observed {T) 
and synthetic {T') spectra at various T^s 



~ \2 



X' 



(2) 



where cr^ is the squared sum of the CALSPEC and our random 
errors, arising primarily from the photometry and to minor extent 
[Fe/H] and log^. Angular diameters are needed to scale syn- 
thetic spectra into physical units: typical 1% internal accuracy 
in fliRFM implies 2% errors in the derived flux. We decided to 
use random errors only because the purpose of the test is exactly 
to verify the range of values allowed once the zero point of the 
temperature scale is assumed. 

Also, the tuning of the absolute calibration in the infrared 
affects the final T^a but it does not modify in any manner the 
shap e of the sy nthe tic spe ctrum, which entirely depends on the 
ICastelli & Kurucj (l2004l) grid interpolated at the proper Teff, 
log^ and [Fe/H]. Notice that we are not searching for the syn- 
thetic spectrum which best matches the observation, rather we 
want to test the effective temperature we derive: while adjust- 
ments to [Fe/H] and \ogg could improve the agreement in the 
blue and visible part, the continuum characteristics are more sen- 
sitive to Teff. 
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Fig. 10. Same as in figure|9]for BD +17 4708. The synthetic spectra have been reddened by E(B-V) - 0.01. Different symbols in the 
right panel correspond to cut longward of 0.50 yum (asterisks) 0.66 yum (diamonds), 0.82 /im (squares) and 1.46 yum (triangles). The 
maximum wavelength used for computing the reduced has been 2 jum to avoid possible contribution from the cool companion. 



4.2.1. HD209458 

For this target we adopted the spectroscopic [Fe/H] - 0.03+0.02 
and \ogg - 4.50 + 0.04 measured fro m the high precision 
HARPS GTO sample (ISousa et al.ll2008l) and used Tycho2 and 
2MASS photometry. We obtain Jeff = 6113 + 49 K, TboI = 
(2.335±0.025)xl0"'^ergcm~2s-' andfl = 0.224 ±0.004 mas in- 
cluding both random and systematic eiTors. The latter result is in 
good agreement with the angular diameters 0.21 5 ± 0.009 mas 
obtained using the new Hipparcos parallaxes (Ivan LeeuwenI 
[2007 ) to convert the linear r adius measured fro m exoplanet 
transit photometry with HST (iBrown et all 1200 ih . Notice that 
~ lOOK cooler effective temperatures would imply values of 
e larger by -3.5% in the IRFM. 

The comparison between the observed and synthetic spec- 
tra at two different Teff is shown in Figure |9l while they both 
succeed to capture the main observed features, the continuum of 
the cooler model is clearly off from the observation. We quan- 
tify the agreement between the HST spectrophotometry and the 
models at various Tgff applying statistics longward of the Ho- 
line (0.66 jjm), the Paschen (0.82 yum) and the Brackett (1.46 yum) 
discontinuity. These cuts define the beginning of the continuum 
in a somewhat arbitrary manner, but they all return consistent 
results thus ensuring that our conclusion is not affected by their 
choice. The reduced is lower than 1 in a roughly ±40 K inter- 
val effectively centered on our prefeiTed solution. While reduced 

< i tells that the size of the errors is still too large to clearly 
favour a solution within that range, the large number of points 
used in the test sets low Icr and 3cr levels, clearly ruling out 
solutions different by ±100 K. 

4.2.2. BD +174708 

This star is the only subdwarf with well measured absolute flux, 
thus making it an important benchmark for testing the tempera- 
ture scale in the metal-poor regime. We adopt the spectroscopic 
parameters [Fe/ H1 = -1.74 ± .09, [a/Fe] = 0.4 and logg = 
3.87 ± 0.08 from lRamfrez et aLl (2006) who also derived Tetf - 
6141 ± 50 K, TboI = (4.89 ± 0.10) x lO^'^ergcm-^s-i and = 
0.1016±0.0023mas. Wecorrectedforreddening£'(B-y) = 0.01 
the optical (TablelTJ and infrared (2MASS) magnitudes, obtain- 
ing Teff = 6120 ± 112 K, Tbo/ = (4.80 ± 0.04) X 10-'* erg cm-2s-i 
and - 0.101 ± 0.003 all in excellent agreement with the afore- 
mentioned analysis. Radial velocities show modulation consis- 
tent with the presence of a low mass companion which could 



influence infrared photometry (iLatham et al.ll98 8^. The flags as- 
sociated with 2MASS indicate excellent quality and no artifact 
nor contamination in any band, pointing toward a negligible ef- 
fect, if any. Nonetheless, since the percent contribution of a cool 
companion increases with increasing wavelength, as safety rule 
we decided not to use Ks in the IRFM though it would change 
the resulting Teff by only 12 K. For our preferred Teff = 6120K, 
shortward of 2 yum there is an outstanding agreement with the 
CALSPEC observed spectrum, meaning that the solution found 
represents well the observation at all wavelengths. A moderate 
increase in the observed with respect to the synthetic flux seems 
to appear longward of 2fim, which could be the signature of 
the cooler companion. On the contrary, cooler solutions overes- 
timate the flux throughout the entire continuum. 

Because of the metal-poor nature of this star, the continuum 
shows up already at bluer wavelengths. We compute the reduced 
X~ in different intervals, starting longward of 0.50 yum: as for the 
previous star, our solution substantially coiTespond to the min- 
ima of all parabolae, independently of the cut adopted. The ran- 
dom errors associated with this star are larger than in the case of 
HD209458, giving shallower minima and thus making it more 
difficult to discriminate between different solutions. However, 
differences up to ±100K are clearly disfavoured (Figure [TOt. 

Summarizing, CALSPEC data support our temperature scale 
which provide the best match to the observed spectrophotometry, 
in both metal-rich and -poor regimes. While differences larger 
than ±40 K are ruled out for HD209458, the observational errors 
for the metal-poor star allow bigger uncertainties. Nonetheless, 
we have determined the fundamental parameters of both stars 
with the same procedure and in both cases our solutions are lo- 
cated at the minimum we regard such a result as a further 
indication that our Tss scale is well calibrated over a wide metal- 
licity range. 

5. The new effective temperature scale 

Our results should be compared with effective temperatures de- 
termined employing different methods. First, we focus on large 
studies which have targeted solar neighbourhood stars, where 
the vast number of objects imposes the use of fast and efficient 
techniques, relying on fitting the observed photometry or spec- 
tra to their synthetic counterpart. An extensive comparison be- 
tween the effective temperatures determined from high resolu- 
tion spectroscopy of solar neighbourhood stars and a version of 
the IRFM similar to that adopted here has been already carried 
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Fig. 11. Upper (lower) panels: comparison between the efFec- 
tiv e temperatures d etermined in this work and those obtained 
bv lValenti & Fische r (2005) (Masana et al. 2006).Areif are this 
- other works in all panels. 



out in ISousa et al ] (I2008h . For metal-poor stars we restrict the 
comparison to purely spectroscopic effective temperatures; their 
validation will be crucial for ongoing and future studies of halo 
stars which are strongly affected by reddening and often lacking 
photometry. 

5.1. Solar Neighbourhood stars 

5.1 .1 . Valenti & Fischer sample 

IValenti & Fischej (12005 ^ have presented a uniform catalogue of 
stellar properties for 1040 nearby F,G and K stars which have 
been observed by the Keck, Lick and AAT planet search pro- 
grams. Fitting the observed spectra with synthetic ones, they 
have obtained effective temperatures, surface gravities and abun- 
dances for every star. For 84 objects in common, there is no ob- 
vious dependence as a function of Teff , except for a drift appear- 
ing below 5000 K. However, when AT^ff is plotted as function 
of metallicity the trend becomes clear, with very significant dis- 
crepancies at the lowest metalhcities (Figure fTTli. 

5.1 .2. Masana et al. sample 

iJVlasana et al.l (l2006h have derived stellar effective temperatures 
and bolometric corrections by fitting V and 2MASS IR photom- 
etry. They calibrate their scale by requiring a set of 50 solar 
analogs drawn from ICavrel de Strobell d 19961) to have on aver- 
age the same temperature as the Sun. 

We have 176 stars in common: there is no obvious trend with 
effective temperatures, and for metalhcities around solar there is 
an overall good agreement. This is not entirely unexpected con- 
sidering that both studies have been calibrated to the Sun (though 
with different approaches): considering [Fe/H] > -1 the mean 
difference (IRFM - Masana) is A^eff = -21 + 6K (cr = 71 K). 
However, when focusing on metal-poor stars [Fe /H] < - 1 there 
is a significantly increasing scatter and a trend resulting in our 
Jeff being cooler up to ~ 200 K at the lowest metalhcities and 
with a mean difference of -95 + 22 K (cr = 157 K). 



5.2. Metal-poor, halo stars 

5.2.1. Temperatures from fits to hydrogen line profiles 

The wings of hydrogen lines are strongly sensitive to the ef- 
fective temperature of the star and only mildly dependent on 
the other stellar parameters, other than being unaffected by red- 
dening. Such approach is particularly effective with metal-poor 
stars, given the lack of severe line blending affecting the hy- 
drogen lines. Thus, provided a proper continuum normaliza- 
tion is applied, wh ich can be non-trivial in some cases (e.g. 
Barklem et al these lines can be used to determine T^g. 

Although significant progress has been made in the last few 
years, the modeling of hydrogen lines (e.g., the Balmer line 
profil es) is still quite uncertain (Barklem et al. 2000; Barklei^ 
I2OO7!) . Nonetheless, the rela tive Teff values der ived in this man- 
ner can be very precise (e.g. lNissen et al.ir2007h . 

We remark that there is no such thing as one Balmer line Teft 
scale, but instead each study depends upon the adopted prescrip- 
tions: LTE vs. NLTE, broadening recipes, mixing-length param- 
eter and even the details on how lines are fitted. Also, the thermal 
structure of the model atmosphere is crucial for the Balmer tem- 
peratures: as concerns ID models, studies relying on OS- instea d 
of ODF-model atmosphere determine hotter T^g ( lGruppl2004l) . 

Aware of the complexity of the picture, in the upper pan- 
els of Figure [12] our IRFM effective temperatures are compared 
with those derived from fits to the Balmer lines in two differ- 
ent studies, which we regard as representative of the LTE and 
NLT E approach, respectiv ely. Circles refer to the comparison 
with iFabbian et al.l (l2009 h who used the Hy6 lines. There is an 
obvious ofi'set, the IRFM returning Teft hotter by 84 + 13 K 
(cr - 66 K), but the small scatter between these two sets fur- 
ther strengthen the conclusion that both techniques have high in- 
ternal precision. A similar conclusion holds also fr om the com- 

S iariso n with the effective temperatures reported in iBergemajuil 
20081 ,and references therein) who used both Ha and HyS fine 
profiles. In this case the difference (IRFM - H lines) is 21 +23 K 
(cr = 72 K) with a possible trend suggesting excellent agreement 
roughly below 6000 K (one star, HD25329 with T^ff = 4785 K 
and AFetf = -15 K is not shown in the upper left panel of Figure 



5.2.2. Excitation equilibrium temperatures 

An important number of iron lines are present in the spectra of 
cool dwarfs, even the metal-poor ones. In an ideal case, the iron 
abundances determined from each of those lines should be con- 
sistent with each other. In practice, however, given an initial set 
of stellar parameters, the line-by-line abundances show trends 
with excitation potential (EP) and/or reduced equivalent width. 
By tuning the stellar parameters, these trends can be eliminated. 
The EP trend is particularly sensitive to Teff , given the strong de- 
pendence of the atomic level populations on temperature, and 
therefore Teff determined by removing the abundance vs. EP 
trend are often refeiTed to as "excitation equilibrium" tempera- 
tures. Because of its nature, this method of Teff determination is 
highly model-dependent. Not only it does require realistic model 
atmospheres and spectrum synthesis, but also accurate atomic 
data and, ideally, a non-LTE treatment of the line formation. The 
advantage of such method is that it is independent of interstellar 
reddening and can be applied to stars with uncertain or unavail- 
able photometry. 

Recently, iHosford et al] (l20()9b have determined LTE excita- 
tion equilibrium temperatures for a sample of metal-poor stars. 
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Fig. 12. Upper panels: comparison between the effective tem- 
pera tures determined in this work and those obtained f rom the 
UB (iFabbian et al.ll2009l circles) and Ha plus li/3 ( Ber gemannI 
Horn squares) line profiles. Lower panel: comparison with re- 
spect to the excitation equilibrium temperatures determined by 



Hosf ord et al. (2009). Two sets of data points are shown because 
Hosford et al., (i2009i) temperatures are sensitive to the uncertain 



logg values of metal-poor stars; squares (triangles) represents 
Teff derived assuming the star to be on the main-sequence (sub- 
giant branch). ATgg are this - other works in all panels. 



The difference found between their temperatures and ours is 
illustrated in Figure [H] (HD 140283 with E{B - V) = 0.000, 
^eff = 5777 K and AT^ff = 8 K is not shown in the lower left 
panel). Because the excitation temperatures are somewhat sensi- 
tive to log g and surface gravities of metal-poor stars are difficult 
to determine due to uncertain/unavailable parallaxes, they pro- 
vide two sets of Teff values, one assuming the star to be on the 
main-sequence (MS) and another one assuming the star to be on 
the subgiant branch (SGB). We remark that for HD 140283 par- 
allax and Balmer jump rule out the main-s equence stage ; our fit 
(Mike Bessell) of the MILES fluxes using iMunari et alj (l2005h 
spectral library provide Tss - 5812/5875Kandlogg - 3.75 for 
E(B - y) = 0.000/0.017, respectively. 

The IRFM temperatures are significantly hotter than the ex- 
citation temperatures by 177 ± 33 K (cr = 122 K) (for their MS 
temperatures) and 240 ± 32 K (cr = 1 16 K) (SGB). In particular, 
the large scatter suggests a decreased relative precision when ap- 
plying excitation equilibrium to very metal-poor stars, so that the 
further investigation of non-LTE effects will be highly desirable 
(Hosford et al. in prep.). 



5.3. The most metal-poor stars in the Galaxy 

Despite theoretical uncertainties on the exact mass range under 
which the first stars formed, it is likely that the most metal-poor 
objects currently observed in the Milky Way halo are second 
generation stars. In case of dwarfs/subgiants, their abundance 
patterns carry dire ct information on the first stars ever formed in 
the Galaxy (e.g. Frebel et al."'2005) and/or on still poorly known 
long time-scale processes which might t ake place below the sur- 
face or deep int o stellar interior (e.g. I Venn & Lambert! l2008t 
iKornet al.ll2009l) . 



Determining their effective temperature and evolutionary 
status (i.e. logg) is crucial to derive reliable abundances and 
constrain different scenarios. At the same time, such a quest 
is in stark contrast with the many practical limitations associ- 
ated with hyper-metal-poor stars: parallaxes are not available 
to help constrain their surface gravities and even when spectra 
with sufficient resolution and S/N are obtained, the model atmo- 
spheres used for the analysis are not yet fully tested at such low 
metallicitie s . Rigor ous analyses s hould also take into account 3D 
dFrebel et al.ll2008h and NLTE ( Aoki et al.ll2006 ') effects, which 
are expected to be considerable in this regime. Determining Teif 
in a way mostly unaffected by the above limitations is not only 
desirable, but also necessary to put spectroscopic analyses on 
firmer grounds. 



5.3.1. HE1 327-2326 

For this star the IRFM returns Teff = 6250 + 60 K in agreement 
within the errors with the spectroscopic value of 61 20 ± 150 K 
obtained from the NLTE analysis of the Balmer lines dKorn et al.l 
|2009), roughly with an offset of the same order of that dis- 
cussed in Section 15.2.11 As we already pointed out, the IRFM 
depends only weakly on the adopted surface gravity: chang- 
ing it by +0.5 dex affects Teff by approximately +2 5 K. In our 
case, w e used logg = 3.7 as recently determined bv iKorn et al.l 
(12009!). The exact metallicity of HE1327-2326 is also uncertain: 
although it is well established that its [Fe/H] < -5.0, estimates 
range from -5.9 to -5.4 depending on the adopted stellar pa- 
rameters and 1D/3D LTE /NLTE analysis performed dAoki et al.l 
120061; iFrebel etalj|2008h . The IRFM is known to depend very 
little on the metallicity and we verified this being particularly 
true (at least in this T^ff regime) for the featureless spectra of 
this hyper-metal-poor star: increasing [Fe/H] by 1 dex in the 
IRFM affects the derived Tss by less then lOK. This conclusion 
supports the suggestion that for hyper-metal-poor stars colour- 
temperature calibration of normal very-metal-poor stars can be 
used instead (see discussion in Section|6]l. 

When running the IRFM for this star we used the n ew grid 
of MARCS model atmosphere (iGustafsson et alJ l2008l) which 
extend down to [Fe/H] = -5.0 and this value was used in 
our implementation. Because of the weak metallicity depen- 
dence discussed above, v ery similar results are obtained if the 
ICastelli & Kuruczld2004 grid (which stops to [M/H] = -4.0) is 
used instead. For the sake of ensuring our results do not depend 
too much on the adopted spectra library, we also checked that 
for stars with higher metallicities MARCS or ATLAS9 models 
return very similar results, with differe nces usually well withi n 
10 K and at most of order 20 K (see also lCasagrande et al.l2006l) . 

We feel the major source of possible systematic error stems 
from reddening, which is very high for this star. We used E{B - 
V) - 0.076 b ased on both extinction maps and in terstellar ab- 
sorption lines dAoki et al.ll2006l: iBeers et al]|200 7') but it should 
be kept in mind that a change of +0.01 mag. in E(B - V) affects 
T^ff by +50 K. 



5.3.2. HE0233-0343 

Though the exact metallicity of this star i s still uncertain, it 
seems well secured as having [Fe/H] < -4.0 dGarcfa Perez et all 
2008, Garcia Perez private communication). Its evolutionary sta- 
tus is also ambiguous, with spectroscopic estimates of logg 
varying from 3.5 to 4.5. Also in this case, the exact values of 
logg and [Fe/H] are not crucial for the IRFM and we checked 




Fig. 14. Upper panels: empirical colour-temperature-metallicity calibrations in the metallicity bins -0.5 < [Fe/H] < 0.5 (filled 
diamonds), -1.5 < [Fe/H] < -0.5 (upward triangles), -2.5 < [Fe/H] < -1.5 (downward triangles) and [Fe/H] < -2.5 (open 
circles). Open squares are for the hyper metal-poor stars HE0233-0343 and HEl 327-2326. Lower panels: residual of the fit as 
function of metallicity. For the two hyper-metal-poor stars, the residual is with respect to the fit at [Fe/H] = -3.5. 
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Fig. 13. Upper left panel: metallicities and effective tempera- 
tures of our sample. All stars have 2MASS and Johnson-Cousins 
photometry. Upper right panel: effective temperatures and gravi- 
ties of our sample. Symbols for different metallicity bins are the 
same as in the left panel. Overplotted fo r reference is a 3 Gyr so- 
lar isochrone from lBertelli et alJ (l2008h . Lower panel: metallic- 
ity sensitivity of our colour-temperature calibration in different 
bands for stars having Teff - 6500 K (top), 6200 K (middle) and 
5900 K (lower) at [Fe/H] = -3.0. 

that changing them even considerably affects T^ff by an amount 
similar to that discussed for HE1327-2326. We adopt [Fe/H] - 



-4.0 and logg - 4.0 from which we derive Teff = 6270 + 
80 K, without acc ounting for possib le systematics arising from 
E(B - y) = 0.025 (iBeers et al.ll2007l) . As we point out in Section 
|6] there might be some issue with the Rc photometry for this 
star. Were we to exclude this band when running the IRFM, T^^ff 
would increase by 25 - 35 K depending on the surface gravity 
assumed. Spectroscopic Teff estimates for this star are still uncer- 
tain, primarily because of its uncertain log^. Were its subgiant 
status to be confirmed, our effective tempe rature would be in 
good agreement with the spectroscopic one (iGarcia Perez et all 
20081) . 

6. Empirical calibrations 

The effective temperatures and the bolometric luminosities de- 
rived via IRFM for our sample allow us to build calibrations 
relating those quantities to the measured colours and metallici- 
ties. As discussed in Section |2] to correctly account for redden- 
ing is crucial though fortunately, for the sake of deriving colour 
relations, reddening affects both the observed photometry and 
the derived fundamental stellar parameters, thus making such 
relations -built using dereddened colours- independent on the 
adopted E(B - V) in first approximation. 

In the following we give the functional form of these cal- 
ibrations, together with the number of stars used, the standard 
deviation obtained in the fitting process and the range of applica- 
bility. The results presented here usually match Casagrande et aT| 
(l2006h within the limits of those calibrations, but extend over a 
wider range now and thus supersede the previous work. Though 
our sample has been assembled explicitly to cover a parameter 
space as large as possible in effective temperature and metallic- 
ity, the detection and observation of stars with [Fe/H] < -2.5 
is still strongly biased around Teff ~ 6500 K. Even if the for- 
mal range of applicability of the calibrations extend well below 
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[Fe/H] < -3, the number of known metal-poor stars consider- 
ably decreases as one moves away from the aforementioned Teff 
(see Figure [TST l. In particular, for metallicities below -4, only 
two stars are currently known, a number clearly inadequate to 
give fits. Fortunately, at these temperatures calibrations at about 
-3.5 seem adequate for even more metal-poor stars, as we dis- 
cuss further in Section 16.11 and 16.21 Nonetheless, we advocate 
particular caution when using these calibrations in poorly sam- 
pled regions of Figure [13] On the contrary for [Fe/H] > -2, 
typical for most of the stellar population observed in the solar 
neighbourhood and Galactic star clusters, our calibrations are 
robust and can be readily used for a number of purposes. 

The core of the present work is to accurately define the zero 
point of the temperature scale in many standard photometric sys- 
tems; we caution however that in some cases real systems might 
not exactly reproduce standard syste ms, especially in the case 
of the faintest sources (Bessell 2005). Users of our calibrations 
should always keep this in mind: although the zero point of the 
Jeff scale is now well defined, in gathering photometry from het- 
erogeneous sources there might be small zero point issues be- 
tween different authors, and this observational uncertainty -if 
present- will introduce small systematic errors to our accurate 
empirical calibrations. 

6.1. Colour-Temperature-Metallicity 

To reproduce the observed Teft versus colour relation and 
take into account the effects of metallicity, the usual fit- 
ting formula has be en adopt ed (e.g. [A lonso et al. 1 996b ; 



Ramirez & Melendez" '2005^ ICasagrande et al.l 120061: 



Gonzalez Hernandez & Bonifacioll2009 ) 



^eff ^ao + aiX + flz^ + a3^[Fe/H] + a4[Fe/H] + fl5[Fe/H]^(3) 

where 0eff - 5040/retf, X represents the colour and a, (/ - 
0, . . . , 5) are the coefficients of the fit obtained iteratively, dis- 
carding points departing more than 3cr. 

The IRFM depends only very mildly on the adopted \ogg 
(Section [2.2b but certain colours could be more affected: for all 
indices we have checked the residual of our calibration and did 
not find any obvious trend with \ogg. Nevertheless, a depen- 
dence on the gravity could be built into the calibrations, since 
log g decreases as one moves from cool dwarfs to hotter turn-off 
stars (Figure [13]). 

The coefficients for various colour indices are given with 
their range of applicability in Table [4] and a comparison between 
the polynomial fits and our sample of stars is shown in Figure 
[T4l We remark that the functional form of Eq. ([3]i may return 
non-physical values when ex trapolated to very low metallicities, 
as exte nsively discussed bvlRvanetalJ (Il999h for the calibra- 
tion of lAlonsoetaD (Il996bi) below FFe/Hl 2.5. We have 

considerably increased the number of very metal-poor (turnoff) 
stars and our calibration behaves as one would expect, i.e. it 
shows a decreasing sensitivity on [Fe/H] when moving from 
-2 to -3, where the metallicity sensitivity vanishes in all bands 
(Figure [T3ll. Moving to [Fe/H] - -4 (or lower), the diverging 
behaviour in Figure [13] reflects the form of the fitting function 
and the values of the coefficients rather than the characteristics 
of metal-poor turnoff stars. In Figure [14] the two hyper metal- 
poor stars (represented by open squares) clearly follow the same 
trend of other iron deficient stars with similar effective temper- 
atures. Using Eq. ([3]l at a fixed [Fe/H] = -3.5 recovers their 
IRFM Teff within the typical accuracy of the calibration. This 
is always true for HEl 327-2326, and also for HE0233-0343 ex- 
cept when using the Rc index, possibly indicating a photometric 



issue in this band for the latter star. This comparison thus war- 
rants the applicability of our calibrations for hyper-metal-poor 
stars if [Fe/H] - -3.5 is assumed and a typical T^ff ~ 6200 K is 
obtained. How well this holds at other effective temperatures is 
still unknown. 

The calibration presented here applies till late K-type 
dwarfs. Those interested in M dwarfs, can instead refer to 
Casagrande et al. (2008): though in that work the zero point has 
not been constrained using solar twins, the absolute calibration 
adopted was similar to that used here, resulting in effective tem- 
peratures approximately on the same scale. Nonetheless, if a link 
between the two scales is needed, we advise users to a careful 
case-by-case study, also considering that the calibration for M 
dwarfs has a different functional form and does not include any 
metallicity term. 

6.1.1. Strbmgren calibration 

The Stromgren index b - y deserves a separate discussion. It is 
often used as a T^ff indicator, but because of its very nature has a 
strong sensiti vity on the m e tallicity and a proper functional form 
is not trivial. lAlonso et al.l (Il996bh excluded the coolest dwarfs, 
where the dependence of b - y upon Teff possibly flattens out. 
Yet, for the most metal poor stars that calibration diverges to 
unphysical values, as discussed in Ryan et al. ( 1999). 

For b - y we have verified that a calibration of the form 
of Eq. ([3]) has strong residuals as function of both colour and 
metallicity and used polynomial fits to correct such trends, i.e. 
Teff = 5O4O/0eff + P{[Pe/li],b - y). To this purpose, we have 
increased the sample with more than 1000 stars from the GCS 
catalogue (iNordstrom et alJ2004l) all having Stromgren photom- 
etry, spectroscopic metallicities from an updated version of the 
Cayrel catalogue (Melendez, in prep.) and for which the IRFM 
could be applied directly using Tycho2 and 2MASS (Casagrande 
et al. in prep.). 

We checked that a third order polynomial in both colour and 
metallicity was enough; the calibration before and after adopt- 
ing such a correction is shown in Figure [15] and the coefficients, 
given in the form P([Fe/H], b-y) = ZU [Fe/H]' +2^^,, Ci(b- 
yY are Mq = -1.9, Mi = 130.4, M2 = 125.7, M3 = 27.4, 
Co = -1003.7, Ci = 7325.9, C2 = -17207.4, C3 = 12977.7. 
Notice that the form of these corrections can lead to unphysical 
values if extrapolated and should never be applied outside of the 
colour and [Fe/H] ranges of Figure [TS] 

6.2. Colour-Flux-Metallicity 

We adopt the same definition of ICasagrande et aTl (l2006l) to de- 
fine the bolometric correction in a given f band, where 



BC( - otboI - nir 



(4) 



and the zero point of the mboi scale is fixed by choosing MboI,© = 
4.74. Empirical bolometric corrections in various bands can thus 
be readily computed using Eq. (|4]i and dereddening the observed 
magnitudes given in Table 8. 

A complementary way of deriving stellar i ntegrated flux via 
photometric indices is given in the form of [Casagrande et al] 
(2006) 

■rBoi(Earth) = lO^^'^'^f^bo + bjX + bzX^ + h^X^ 

+ b4X[Fe/n] + bslFe/U] + bdFs/Uf^. (5) 
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Table 4. Coefficients and range of applicability of the colour-temperature-metallicity relations. The photometric systems are 
Johnson-Cousins BV{RI)c, 2MASS JHKs, Tycho2 {BV)t and Stromgren by. For the latter, additional coiTections as function 
of [Fe/H] and (b - y) apply, as discussed in Section 16.1.11 For some indices the calibrations are given down to [Fe/H] - -5.0, 
meaning that the effective temperatures of such a metal-poor star can be recovered using [Fe/H] = -3.5 in Eq. Notice that only 
two hyper metal-poor stars are currently known and caution should be used, as discussed in the text. Especially for metal-poor stars, 
please refer to Figure[T3]to check that the calibration is not extrapolated outside its [Fe/H] range. 
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N is the number of stars employed for the fit after the 3cr clipping and cr(T^ff) is the standard deviation (in Kelvin) of the proposed calibrations. 
Notice that the standard deviation does not account for the uncertainty in the zero point of the temperature scale, which is of order 15 - 20 K 
(Section[32]and|4TJ. 




As for the temperature calibrations, also in this case the 
fluxes of the two hyper metal-poor stars can be recovered adopt- 
ing [Fe/H] = -3.5 in Eq. (|5]l, though we caution that the license 
of this approach for considerably bluer or redder indices is still 
unknown. 



6.3. Colour-Angular diameters 

Limb-darkened angular diameters can be readily derived from 
the basic definition involving effective temperatures and bolo- 



metric fluxes, using the calibrations given in Sections 16.11 
and 16.21 Nonetheless, very tight and simple relations exist in 
the J band and in Ta ble |6] we give them in the form of 
ICasagrande et al.1 (l2006h 

e^CQ + ci -^(p{mj,X) (6) 
where 

0(my,X) = lO""'*"'^^ (7) 

for a given colour index X. These relations show remarkably 
small scatter and no metallicity dependence, thus proving ideal 
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Fig. 15. Upper panel: empirical colour-temperature-metallicity 
calibration in b - y before (dotted) and after (continuous lines) 
the polynomial correction. Central and lower panels: residuals 
before and after the polynomial corrections. 



to build a network of small calibrators for interferometric mea- 
surements, for characterizing extrasolar planet transits or mi- 
crolensing events. 

6.4. The colours of the Sun 

The interpolation of Eq. CJ at Teft = 5777 K and [Fe/H] = 
returns the colours of the Sun, which are given in Table Q For 
the Tycho2 and 2MASS system, those can be readily compared 
with the averaged ones from the twins of Section 13.11 not un- 
expectedly there is good agreement, all but one within few mil- 
limag, which usually (at maximum) correspond to few (20) K in 
Jeff- We have also checked that fitting our twins as function of 
Jeff and [Fe/H] returns colours almost identical to their average, 
further confirming that our sample of twins is homogeneously 
distributed in temperature and metallicity around the colours of 
the Sun inferred from our scale. 

In recent years, there has been cons iderable work in order 



to determine the colours of th e Sun (e . g. IS ekiguch i & Fukueital 
2000; Ramirez & Me lendezI l2005bl: iHolmberg et alj l2006t 
PasQuini et al. 2008). One of the most extensive analysis is that 
of Holmberg et al. (2006): the remarkably good agreement we 
have in the optical colours can be understood from the depen- 
dence of these indices on both T^ff and [Fe/H]. The approx- 



imately 1 00 K cooler effective temperature scale adopted by 
IHolmberg etlD (l2006h favours bluer colours, which are grossly 
compensated to the red by the underestimation of ~ 0.1 dex 
in the GCS photometric metallicities with respect to spectro- 
scopic ones selected to be consistent with our temperature scale 
(IHolmberg et al.ll2009l) . Our B -V ^ 0.641 is also in very good 
agreement with the B - V - 0.649 + 0.016 found studying so- 
lar twins in M67 (Pasquini et al. 2008). For this cluster, using 
our colour-temperature relation to compare V - Ks photometry 
with theoretical isochrones shows remarkably good agreement 
(Vandenberg, private communication). 

Infrared indices derived inverting Eq. (O depend almost ex- 
clusively on the adopted Tgs scale, which is respon s ible fo r our 
much redder colours than those of IHolmberg et alj (|2006|) . Our 
V- J, H and Ks are in good agreement with those reported in 
i Rieke et al. (2008) and obtained from solar- type stars or com- 
puted convolving various solar spectra with the appropriate filter 
curves and using their revised absolute physical calibration. 

The empirical colours in Table [T] are also in agreement with 
the synthetic ones, computed using the same zero points and ab- 
solute calibration for Vega used in the IRFM to derive our T^s 
scale. Therefore, the uncertainty in the zero points used to gen- 
erate synthetic colours is at the smallest level possible, yet of the 
order of 0.01 mag. (Section [3 .2l l. allowing us to address the re- 
liability of the models at this level of precision. While using a 
theoretical spectra of Vega may (partly) compensate model in- 
accuracies in the process of setting the zero points, the approach 
adopted here allows us to focus on the quality of the solar syn- 
thetic spectra. The agreement is remarkable, on the order of 0.01 
mag. and never exceeding 0.02, which is also of the same size of 
the difference between those synthetic models. 

7. Conclusions 

The primary goal of this work has been to provide a new absolute 
effective temperature scale. An unprecedented accuracy of few 
tens of Kelvin in the zero point of our scale has been achieved 
using a sample of solar twins. For these stars the high degree of 
resemblance to the Sun has been determined entirely model in- 
dependently, without any prior assumption on their physical pa- 
rameters, most importantly T^g. Notice that by calibrating our 
results via solar twins we are entirely unaffected from possi- 
ble issues and uncertainties related to Vega. Nonetheless, we 
regard as comforting that our findings are in close agreement 
with the lat est absolute fluxes dCohen et al.ll2003l: lBohlinll2007l: 
iRieke et ani2008i) . We further took advantage of such a promis- 
ing situation by fine-tuning the adopted fluxes so as to improve 
the consistency of the effective temperatures determined from 
each band used in the IRFM. This methodology gives us con- 
fidence that the stellar parameters we determined are well cal- 
ibrated not only around the solar value, but over a wide range 
in Teff and [Fe/H]. Notice that the IRFM is little model depen- 
dent and certainly not at the solar value because of our calibra- 
tion procedure. Small spurious trends arising from the adopted 
library at different temperatures and metallicities can not be en- 
tirely ruled out, but should be small. Though the zero point of 
our new Tetj scale is entirely set by solar twins, it agrees within 
few degrees with independent verifications conducted via inter- 
ferometric angular diameters and HST spectrophotometry in the 
metal-poor and -rich regimes. 

In the process of establishing the zero point of the effective 
temperature scale via IRFM, we nailed down the differences 
with respect to other implementations of the same technique. 
We have used two independent IRFM versions to study the dis- 



Casagrande et al.: The effective temperatiire scale: resolving different versions 



19 



crepancies among various temperature scales that appeared in 
literature over the years and proved that the absolute caUbration 
of the photometric systems used was responsible for explain- 
ing most of the differences. At solar temperatures and metallici- 
ties the long-standing dichotomy between photometric and spec- 
troscopic Teff is easily explained once it is understood that the 
IRFM can in principle accommodate any temperature scale since 
its zero point depends on the absolute calibration of the photom- 
etry adopted. The main goal of the present paper has been ex- 
actly to tackle this issue using the best constraint available to 
date. 

The improved bolometric fluxes determined for metal-poor 
stars have also been used to put on firmer ground the temper- 
ature scale in this rather unexplored regime. For metallicities 
typical of halo stars our T^g scale is roughly 100 K hotter than 
those determined from the Balmer lines and 200 K hotter than 
those obtained from the excitation equilibrium. While spectro- 
scopic effective temperature determinations have considerable 
model dependence and are degenerate with other stellar parame- 
ters (namely log g and [Fe/H]), the IRFM offers a powerful alter- 
native, free from any of the above limitations. However, relying 
on the photometry, the IRFM is influenced by reddening, which 
becomes a considerable source of uncertainty when targeting ob- 
jects outside of the local bubble. For our sample of metal-poor 
stars we have been cautious in determining reddening as best 
we could. Our improved determination of E{B - V) also explain 
the remaining discrepancies with other T^ff scales. We think the 
effective temperatures determined for our sample of stars will 
serve to better calibrate spectroscopic T^s determinations. This 
will be particularly relevant when large spectroscopic surveys 
targeting different stellar populations in the Galaxy start operat- 
ing: support from the existing or forthcoming photometric sur- 
veys will be possible only if reddening will be determined on a 
star-by-star basis. We feel this will not be possible in many cases 
and stellar parameters will have to rely on spectroscopy only. 

Based on our sample of dwarfs and subgiants, a set of ho- 
mogeneously calibrated colours versus temperatures, bolomet- 
ric fluxes and angular diameters have also been determined. A 
number of problems of interest to stellar and Galactic Chemical 
evolution depend on the assumption made in these relations and 
our results will permit those problems to be tackled with greater 
confidence. 
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-0.1789 
-0.2609 
-0.2458 
-0.2407 
-0.1276 
-0.3756 
-1.1377 
-0.5392 
-0.4354 
-0.1301 
-0.2929 
-0.3935 
-0.2339 
-0.4491 
-0.2991 
-0.3112 
-0.2147 
-0.2506 
-0.1229 
-0.1968 
-0.3054 
-0.2408 
-0.2508 
-0.2604 
-0.2509 
-0.0358 
-0.1418 
-0.2329 
-0.1866 
-0.1908 
-0.2132 
-0.1891 
0.0132 
-0.0769 
-0.1854 
-0.1232 
-0.1059 
-0.1577 
-0.1326 
0.1274 
-0.0012 
-0.1169 
-0.0512 
-0.0098 
-0.0665 
-0.0407 
0.1098 
-0.0186 
-0.1469 
-0.0690 
-0.0315 
-0.0879 
-0.0641 



-0.0720 
0.0343 
0.0809 
0.0669 
-0.0413 
-0.0048 
0.0019 
-0.0011 
-0.0432 
-0.0067 
0.0411 
0.0089 
0.0047 
0.0346 
0.0240 
-0.0420 
-0.0252 
-0.0166 
-0.0231 
-0.0213 
-0.0124 
-0.0160 
-0.0142 
-0.0186 
-0.0171 
-0.0187 
-0.0181 
-0.0185 
-0.0184 
-0.0015 
-0.0112 
-0.0123 
-0.0122 
-0.0118 
-0.0138 
-0.0123 
0.0057 
-0.0048 
-0.0095 
-0.0066 
-0.0054 
-0.0092 
-0.0078 
0.0179 
0.0022 
-0.0047 
-0.0009 
0.0031 
-0.0016 
0.0005 
0.0163 
0.0006 
-0.0066 
-0.0019 
0.0013 
-0.0029 
-0.0017 



260 
261 
255 
262 
253 
249 
261 
256 
331 
186 
202 
196 
332 
328 
363 
241 
177 
197 
184 
314 
369 
316 
179 
180 
203 
185 
184 
196 
193 
200 
191 
177 
197 
186 
187 
193 
346 
186 
195 
192 
308 
303 
314 
331 
184 
192 
195 
353 
344 
363 
353 
190 
199 
201 
328 
317 
328 



3.1 
4.2 
4.0 
3.9 
2.7 
0.7 
0.9 
0.9 
2.3 
1.9 
3.2 
1.6 
2.8 
2.9 
2.8 
1.9 
0.7 
1.5 
0.8 
0.8 
0.9 
0.9 
0.8 
0.7 
1.0 
0.7 
0.7 
0.7 
0.7 
1.0 
1.1 
0.9 
0.8 
1.1 
1.0 
0.9 
3.4 
2.9 
3.0 
2.6 
0.8 
1.9 
1.7 
4.6 
3.6 
3.6 
3.1 
2.8 
1.0 
2.2 
4.9 
3.8 
3.6 
3.2 
2.9 
2.5 
1.1 



N is the number of stars employed for the fit after the 3 sigma clipping and cr(%) is the standard deviation of the final caUbrations in percent. 
The coefficients of the calibrations bi are given in units of lO^^'^ergcm^^ s"'. 
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Table 6. Coefficients and range of appUcability of the angular diameter caUbrations. 



Colour 


Colour range 


Co 


Cl 


N 


o-(%) 


V- 


J 


[0.73,2.44] 


0.00015 


4.65293 


394 


2.2 


V - 


H 


[0.88,3.01] 


0.00004 


4.15613 


389 


1.8 


V - 


■Ks 


[0.93,3.15] 


0.00020 


4.05037 


394 


1.9 


Vt 


-J 


[0.72, 2.56] 


0.00218 


4.44568 


270 


2.1 


Vt 


-H 


[0.87,3.14] 


0.00227 


3.98945 


255 


1.5 


Vt 




[0.92, 3.27] 


0.00286 


3.88433 


268 


1.5 



N is the number of stars employed for the fit after the 3it clipping and it(%) is the standard deviation of the calibrations. 



Table 7. The colours of the Sun. For the indices obtained inverting our T^fi scale, random eiTors are from the dispersion of the fits 
in Table m The uncertainty on the zero point of our T^ff scale is of order 15 K (Section [372b . which usually implies systematic eiTors 
considerably smaller than the random ones. The only exception is for optical-infrared indices which are very sensitive to Teff and 
show small intrinsic scatter, of the order of the aforementioned zero point uncertainty. Also shown for comparis on are the averaged 
colours and standar d deviation of the solar twins, as well as the synthetic colours computed from the ATLAS9 dCastelli & Kuruca 
l2004i) and MARCS dGustafsson et al.ll2008l) models. 







Teff scale 
(rand. + syst. errors) 


MARCS 


ATLAS9 


Twins 


B- 


- V 


0.641 ±0.024 ±0.004 


0.622 


0.645 




V- 


Rc 


0.359 ±0.010 ±0.003 


0.357 


0.358 




(R- 


De- 


0.333 ±0.010 ±0.002 


0.347 


0.349 




V- 


le 


0.690 ±0.016 ±0.004 


0.704 


0.707 




V- 


-J 


1.180 ±0.021 ±0.007 


1.171 


1.180 




V- 


H 


1.460 ±0.023 ±0.010 


1.458 


1.479 




V - 


Ks 


1.544 ±0.018 ±0.010 


1.543 


1.553 




J- 


Ks 


0.362 ± 0.029 ± 0.003 


0.372 


0.373 




(B- 


V)t 


0.730 ±0.031 ±0.006 


0.723 


0.750 


0.735 ± 0.024" 


Vt 


-J 


1.254 ±0.022 ±0.008 


1.240 


1.250 


1.254 ±0.041" 


Vt - 


-H 


1.534 ±0.019 ±0.011 


1.527 


1.550 


1.549 ±0.048" 


Vt - 


-Ks 


1.619 ±0.013 ±0.011 


1.612 


1.623 


1.623 ±0.040" 


b- 


-y 


0.409 ±0.010 ±0.002 






0.409 ±0.003* 



" Average and standard deviation of the colours in Table[2] 

* From Melendez et al. in prep., fitting solar twin colours as a function of T^s and [Fe/H]. 
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Appendix A: Comparing the TCS and 2I\/IASS 
absolute calibration 



The absolute calibration of lAlonso et al. was ob- 

tained applying the IRFM to a sample of stars for which 
direct measurements of angular diameters were available. 
Because of the difficulties involved in achieving milli- 
arcsecond resolution, that sample was almost entirely com- 
posed of giants with angular diameters measured via Lunar 
Occultations and Michelson Interferometry (Teff < 5000 K) 
or Intensity Interferometry (Teff > 6000 K). One of 
the intriguing results of that analysis was the impossibil- 
ity of setting the same zero point of the absolute calibra- 
tion using ang ular diameters m easured b y Lunar Occultations 
dWhite & Feierman.1987: .Ridgwav et al.lll980.) and Michelson 



Table A.l. Characteristic parameters of the 2MASS and TCS 
photometric systems. Wavelengths are in A and the Vega's 
monochromatic absolute fluxes in erg cm"^ s"' A"' . 



Interferometry (iHutter et al.lll989H di Benedetto & Rabbiall987l; 
iMozurkew ich et al.' '1991') with those measured by Intensity 
Interferometry (Hanburv Brown et al. 1974). The absolute cal- 
ibration (in the Johnson system) proposed by Alonso et al.l 
(Il994ah is a weighted average from their table 10 and it 
is interesting to notice that the one derived from Intensity 
Interferometry alone is 4.8 (J) 1.3 (H) and 4.0 (K) percent lower 
than the averaged, proposed one. As we have discuss through- 
out the papeiQ, lower infrared fluxes support higher T^ff (in this 
case, the average difference in Johnson system would be 3.4% 
supporting T^ff approximately hotter by 70 K), so it is not sur- 
prising our effective temperature scale provides good agreement 
with interferometric measurements despite being considerably 
hotter than most of the previous IRFM analyses. 

To gauge a further insigh t into the problem, here we directly 
compare the TCS (given in lAlonso et al. 1994b) and 2MASS 
absolute calibration. Such an exercise, however is not straight- 
forward since the absolute calibration in different photometric 
system is obtained using different filter transmission curves and 
therefore is associated to different effective wavelengths. In ad- 
dition, for the sake of the IRFM, in any given band ^, it is the 
composite effect of Vega's magnitudes and fluxes which mat- 
ters, i.e. 10" '*™f. Therefore, for a meaningful comparison we 
need to refer everything to a common wavelength, the 2MASS 
one being the natural choice in this case. This is done in Table 
lA.ll bv computing ^eff i-e the composite effect of magnitudes 
and fluxes shifted to the 2MASS Asg in the case of TCS (Figure 

The 2MASS absolute calibration is on the average lower than 
the TCS by 4.6% (a value qualitatively in agreement with the dif- 
ference in the Johnson system discussed above), thus returning 
Teff on average hotter by ~ 90 K, and explaining the bulk of the 
differences discussed in Section 12.31 when comparing the sam- 
ple stars directly. Similar conclusions can be drawn when com- 
paring with the absolute fluxes and magnitudes of Vega used in 
table 1 of Gonzalez Hernandez & Bonifacio (2009). In this case 
the photometric system is the same (2MASS) and one can di- 
rectly compare ^Fifj: the difference is -3.3% (/), H-1.3% (H) and 
-2.8% (Ks) thus giving an average of -1.6% which correspond 
to ~ 30 K, again in line with the differences discussed in Section 





2MASS 


TCS 


A(%) 




12285 


12790 




mj 


-0.001 


-0.013 




Ty 


3.079 


2.912 




^ff 


3.076 


3.303 


-6.9% 




16385 


16483 




niH 


+0.019 


-0.005 




Th 


1.150 


1.192 




^ff 


1.170 


1.211 


-3.4% 


■^eS 


21521 


21869 




IHk 


-0.017 


-0.029 




Tk 


0.430 


0.426 




'Fee 


0.423 


0.439 


-3.6% 



2M ASS effective wav e length s {A^s), magnitudes and fluxes are 
from lCasagrande et all ( I2006I) . where the latter has been modified 
by -1.6, +1.5 and +0.3 percent in J,H and Kg ban d respectively 
as desc ribed in Section[3j2] The TCS values are from lAlonso et al.l 
(11994b '). For the TCS system T^s has been computed by shifting 
the value at the 2MASS effective wavelength (Figure lATt . A(%) is 
the percent decrease of the TCS Tss needed to match the 2MASS 
values. 



We have verified using our IRFM implementation that a 1% in- 
crease in infrared fluxes correspond to a decrease of 20 K in T^g. 
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1.10x10'' 1.25x10'' 1.40x10'' l.SQxlO* 1.65x104 1.80x10* 1.90x10'' 2.15x10'' 2.40x10* 

X (A) X (A) X (A) 



Fig. A.l. Comparison between TeS in 2MASS (filled squares) and the TCS before (open circles) and after (filled circles) correcting 
for the same /igif • The correction has been done shifting the TCS values along the continuum of Vega, obtained by fitting a second 
order polynomial to the observed spectral energy distribution (from Bohlin 2007.) . Overplotted for comparison are the 2MASS 
(dotted lines) and TCS (continuous lines) filter transmission curves. 



Table 8. Sample stars". 



name 




[Fe/H] 


mBol 

















B 


V 


Rc 


Ic 


J 


H 


Ks 


E(B - V) 


BD+000444 


4, 


.70 


-0 


.01 


8.802 


4340 


+ 


74 





.254 


± 





.009 


10.732 


9.532 


8.796 


8.149 


7.260 


6.692 


6.500 


0, 


.000 


BD+002058 


4, 


.19 


-1, 


.16 


10.020 


6245 


+ 


58 





.070 


+ 





.001 


10.652 


10.199 


9.890 


9.570 


9.165 


8.899 


8.842 


0, 


.017 


BD+004470 


4. 


.61 


-1 


.53 


10.351 


6358 


+ 


61 





.058 


+ 





.001 


11.040 


10.450 








9.495 


9.266 


9.207 


0, 


.000 


BD+012831 


4, 


.50 


-1 


.27 


10.677 


5331 


+ 


50 





.071 


+ 





.001 


11.549 


10.882 








9.455 


9.063 


8.993 


0, 


.000 


BD+013597 


4, 


.04 


-1, 


.87 


10.809 


6289 


+ 


63 


0, 


.048 


+ 


0, 


.001 


11.544 


11.089 








10.005 


9.718 


9.665 


0, 


.045 


BD+023375 


4, 


.12 


-2, 


.24 


9.699 


6102 


± 


88 





.085 


± 





.003 


10.414 


9.944 


9.635 


9.297 


8.799 


8.537 


8.498 


0, 


.022 


BD+024651 


3, 


.79 


-1 


.78 


9.980 


6349 


+ 


67 





.069 


+ 





.002 


10.646 


10.204 








9.193 


8.902 


8.831 


0, 


.038 


BD+030443 


4, 


.50 


-0 


.93 


9.417 


4591 


+ 


66 





.171 


+ 





.005 


10.890 


9.960 








7.940 


7.428 


7.295 


0, 


.000 


BD+030740 


3, 


.97 


-2, 


.70 


9.584 


6419 


+ 


83 





.081 


+ 





.002 


10.168 


9.808 








8.795 


8.532 


8.516 


0, 


.022 


BD+040415 


4, 


.50 


-0, 


.63 


9.435 


4833 


± 


50 


0, 


.153 


± 


0, 


.003 


10.737 


9.825 


9.270 


8.783 


8.045 


7.553 


7.449 


0, 


.000 


BD+053640 


4, 


.56 


-1 


.22 


10.154 


5092 


± 


57 


0, 


.099 


+ 





.002 


11.160 


10.430 








8.854 


8.450 


8.341 


0, 


.000 


BD+054481 


4, 


.50 


-0 


.98 


9.975 


5727 


+ 


65 





.085 


+ 





.002 


10.747 


10.130 








8.884 


8.581 


8.537 


0, 


.000 


BD+072634 


4, 


.50 


-0 


.66 


9.595 


5469 


+ 


64 





.111 


+ 





.003 


10.420 


9.780 








8.407 


8.095 


7.980 


0, 


.000 


BD+074841 


4, 


.09 


-1 


.39 


10.260 


6133 


+ 


61 





.065 


+ 





.001 


10.840 


10.380 








9.334 


9.054 


9.017 


0, 


.000 


BD+083095 


4, 


.13 


-0, 


.78 


9.855 


5690 


± 


58 


0, 


.091 


± 


0, 


.002 


10.586 


9.996 








8.750 


8.410 


8.375 


0, 


.000 


BD+090352 


4, 


.32 


-2, 


.05 


9.970 


6145 


+ 


98 





.074 


± 





.002 


10.616 


10.174 








9.012 


8.829 


8.749 


0, 


.012 


BD+092190 


4, 


.03 


-2, 


.54 


10.971 


6477 


+ 


75 





.042 


+ 





.001 


11.533 


11.147 








10.193 


9.932 


9.907 


0, 


.015 


BD+092879 


4, 


.99 


-0 


.05 


10.239 


4895 


+ 


41 





.103 


+ 





.002 


11.588 


10.600 


10.034 


9.563 


8.864 


8.368 


8.272 


0, 


.000 


BD+095076 


4, 


.50 


-0 


.08 


9.340 


5360 


+ 


49 





.130 


+ 





.002 


10.280 


9.514 








8.118 


7.748 


7.627 


0, 


.000 


BD+ 100449 


4, 


.40 


-1, 


.75 


10.581 


5832 


± 


60 


0, 


.062 


± 


0, 


.001 


11.265 


10.757 


10.417 


10.067 


9.573 


9.287 


9.180 


0, 


.000 


BD+1 10299 


4, 


.50 


-0 


.87 


10.474 


5929 


± 


55 





.063 


+ 





.001 


11.160 


10.600 








9.491 


9.174 


9.124 


0, 


.000 


BD+1 10468 


4, 


.58 


-1 


.55 


10.582 


5739 


+ 


50 





.064 


+ 





.001 


11.312 


10.774 








9.549 


9.229 


9.135 


0, 


.004 


BD+1 12021 


4, 


.50 


0.09 


9.453 


5262 


+ 


56 


0, 


.128 


+ 





.003 


10.414 


9.632 


9.192 


8.792 


8.190 


7.805 


7.689 


0, 


.000 


BD+1 14725 


3, 


.59 


-0 


.85 


9.384 


5476 


± 


45 





.122 


± 


0, 


.002 


10.193 


9.548 








8.223 


7.856 


7.757 


0, 


.000 


BD+132567 


4, 


.00 


-1 


.18 


9.642 


5386 


± 


57 





.112 


± 





.002 


10.481 


9.828 


9.438 


9.038 


8.431 


8.051 


8.002 


0, 


.000 


BD+1 32698 


3, 


.98 


-0 


.94 


9.245 


5824 


+ 


54 





.115 


+ 





.002 


9.956 


9.376 








8.187 


7.900 


7.844 


0, 


.000 


BD+133683 


3, 


.93 


-2, 


.38 


10.060 


5488 


+ 


62 





.089 


+ 





.002 


11.217 


10.566 


10.100 


9.640 


9.005 


8.591 


8.527 


0, 


.080 


BD+141947 


4, 


.50 


-0 


.42 


9.952 


5131 


+ 


51 





.107 


+ 


0, 


.002 


10.980 


10.180 








8.659 


8.200 


8.137 


0, 


.000 


BD+150028 


4, 


.39 


-0 


.72 


9.633 


5891 


± 


50 





.094 


± 





.002 


10.302 


9.743 








8.593 


8.313 


8.245 


0, 


.000 


BD+174708 


3, 


.87 


-1 


.74 


9.312 


6120: 




112 





.101 


+ 





.003 


9.922 


9.476 


9.183 


8.854 


8.435 


8.108 


8.075 


0, 


.010 


BD+183423 


4, 


.17 


-0 


.88 


9.679 


6059 


+ 


58 





.087 


+ 





.002 


10.266 


9.783 








8.682 


8.448 


8.384 


0, 


.000 


BD+191730 


4, 


.40 


-2, 


.21 


10.568 


6210 


+ 


65 





.055 


+ 





.001 


11.178 


10.736 


10.436 


10.136 


9.692 


9.414 


9.396 


0, 


.000 


BD+ 192646 


4, 


.50 


-0, 


.75 


9.732 


5640 


± 


58 


0, 


.098 


± 


0, 


.002 


10.527 


9.887 








8.619 


8.268 


8.231 


0, 


.000 


BD+202030 


4, 


.00 


-2, 


.64 


10.958 


6496 


± 


65 


0, 


.042 


± 





.001 


11.580 


11.200 








10.220 


9.980 


9.937 


0, 


.024 


BD+202594 


4, 


.49 


-0 


.87 


9.874 


6041 


+ 


52 





.080 


+ 





.002 


10.479 


9.992 








8.894 


8.642 


8.576 


0, 


.000 


BD+202972 


4, 


.00 


-0 


.36 


10.340 


5681 


+ 


46 





.073 


+ 





.001 


11.070 


10.450 








9.220 


8.890 


8.826 


0, 


.000 


BD+2 12244 


4, 


.50 


-1 


.01 


9.794 


5504 


+ 


58 





.100 


+ 





.002 


10.560 


9.960 








8.665 


8.261 


8.192 


0, 


.000 


BD+224454 


4, 


.49 


-0, 


.50 


9.301 


5326 


± 


62 


0, 


.134 


± 


0, 


.003 


10.270 


9.500 








8.048 


7.707 


7.611 


0, 


.000 


BD+241676 


3, 


.84 


-2, 


.54 


10.610 


6387 


+ 


92 





.051 


+ 





.002 


11.171 


10.814 






9.811 


9.545 


9.541 


0, 


.013 


BD+244460 


4, 


.38 


-0 


.89 


9.223 


5160 


+ 


62 





.148 


+ 





.004 


10.235 


9.463 


9.018 


8.609 


7.931 


7.513 


7.461 


0, 


.000 


BD+251981 


3, 


.87 


-1 


.34 


9.259 


6918 


+ 


88 





.081 


+ 





.002 


9.627 


9.320 






8.592 


8.438 


8.392 


0, 


.000 


BD+262251 


4, 


.18 


-0 


.95 


10.254 


6049 


+ 


46 





.067 


+ 





.001 


10.846 


10.370 






9.290 


9.022 


8.963 


0, 


.000 


BD+262606 


4, 


.18 


-2, 


.36 


9.534 


6194 


± 


56 


0, 


.089 


± 


0, 


.002 


10.157 


9.732 






8.676 


8.394 


8.352 


0, 


.005 


BD+262621 


4, 


.00 


-2, 


.54 


10.822 


6336 


± 


72 





.047 


± 





.001 


11.428 


11.023 






9.994 


9.749 


9.731 


0, 


.009 



n 



a" 



O 

s 

I 

a' 



BD+263578 


3, 


.88 


-2, 


.29 


9.188 


6425: 


t 


104 





.097 


+ 





.003 


9.749 


9.360 
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8.802 


8 


.357 


7.966 


7.354 


6.971 


6.886 


0. 


.000 


O 


HD88725 


4, 


.37 


-0.64 


7.633 


5737 + 


52 





.249 


+ 





.005 


8.358 


7.753 









6.543 


6.242 


6.153 


0. 


.000 




HD89813 


4, 


.51 


-0.17 


7.612 


5379 + 


53 





.286 


+ 





.006 


8.553 


7.769 


7, 


.346 


6.963 


6.399 


6.032 


5.911 


0. 


.000 




HD90508 


4, 


.35 


-0.33 


6.319 


5743 + 


70 


0, 


.455 


+ 


0, 


.012 


7.049 


6.445 









5.195 


4.895 


4.874 


0. 


.000 


& 


HD90663 


4, 


.63 


-0.29 


8.203 


4943 + 


41 





.258 


+ 





.004 


9.450 


8.530 


7 


.989 


7.550 


6.851 


6.386 


6.265 


0. 


.000 




HD92786 


4, 


.64 


-0.29 


7.835 


5309 + 


60 





.265 


+ 





.006 


8.770 


8.020 


7, 


.593 


7.208 


6.572 


6.227 


6.133 


0. 


.000 


B 


HD92788 


4, 


.42 


0.28 


7.229 


5717 + 


64 





.302 


+ 





.007 


8.000 


7.310 









6.131 


5.798 


5.721 


0. 


.000 


< 


HD94028 


4, 


.36 


-1.44 


8.064 


6110 + 


66 





.180 


+ 





.004 


8.640 


8.202 


7, 


.917 


7.585 


7.130 


6.854 


6.832 


0. 


.000 


5' 


HD9562 


4, 


.02 


0.19 


5.704 


5876 + 


71 


0. 


.577 


+ 


0, 


.015 


6.397 


5.757 


5, 


.407 


5.087 


4.627 


4.391 


4.258 


0. 


.000 




HD97320 


4, 


.23 


-1.28 


8.047 


6134 + 


55 





.180 


+ 





.004 


8.642 


8.161 


7, 


.866 


7.554 


7.137 


6.868 


6.790 


0. 


.000 


o 


HD97916 


4, 


.02 


-0.94 


9.126 


6451 + 


69 





.099 


+ 





.002 


9.628 


9.211 




— 





8.319 


8.104 


8.018 


0. 


.000 




HD98281 


4, 


.61 


-0.20 


7.123 


5426 + 


56 


0, 


.352 


+ 


0, 


.008 


8.022 


7.275 


6 


.862 


6.487 


5.914 


5.575 


5.457 


0. 


.000 




HD99109 


4, 


.32 


0.38 


8.912 


5266 + 


46 





.164 


± 





.003 


9.970 


9.100 









7.626 


7.259 


7.162 


0. 


.000 


1 


HD99747 


4, 


.12 


-0.50 


5.799 


6621 + 


76 





.435 


+ 





.011 


6.221 


5.857 









4.997 


4.825 


4.797 


0. 


.000 




HE0024-2523 


4, 


.30 


-2.67 


14.791 


6583 + 


99 





.007 


+ 





.001 


15.321 


14.913 









14.038 


13.748 


13.749 


0. 


.006 




HEO 130-2303 


4, 


.60 


-2.90 


14.497 


6590 + 


134 





.008 


+ 





.001 


15.141 


14.760 


14.505 


14.207 


13.895 


13.550 


13.589 


0. 


.004 




HE013 1-2740 


3, 


.03 


-3.08 


14.412 


5487 + 


61 


0, 


.012 


+ 


0, 


.001 


15.214 


14.630 


14.240 


13.836 


13.246 


12.865 


12.788 


0. 


.013 




HE0148-2611 


4, 


.30 


-3.00 


14.238 


6595 + 


91 





.009 


+ 





.001 


14.824 


14.453 









13.559 


13.298 


13.303 


0. 


.004 




HE02 18-2738 


4, 


.30 


-3.52 


14.792 


6581 + 


81 





.007 


+ 





.001 


15.277 


14.883 









13.968 


13.762 


13.699 


0. 


.004 




HE0233-0343 


4, 


.00 


-4.30 


15.337 


6270 + 


80 





.006 


+ 





.001 


15.867 


15.430 


If 


i.l71 


14.818 


14.359 


14.077 


14.063 


0. 


.025 




HEl 148-0037 


3, 


.70 


-3.46 


13.446 


6345 + 


74 





.014 


+ 





.001 


14.018 


13.614 








12.610 


12.318 


12.272 


0. 


.022 




HE1327-2326 


3, 


.70 


-5.00 


13.090 


6250 + 


60 


0, 


.017 


+ 


0, 


.001 


14.016 


13.535 


13 


i.211 


12.854 


12.357 


12.068 


11.986 


0. 


.076 




HE2133-1426 


4, 


.10 


-2.80 


15.345 


6259 + 


97 





.006 


± 





.001 


15.935 


15.486 


If 


i.l91 


14.845 


14.380 


14.139 


14.111 


0. 


.028 





LP635-14 
LPS 15-43 
LP831-70 



4.03 

4.19 
4.38 



-2.48 

-2.80 
-2.93 



10.994 

10.688 
11.413 



6520 + 70 

6535 + 68 
6408 + 67 



0.041 + 0.001 11.760 11.330 11.050 10.740 10.285 10.041 9.997 

0.047 + 0.001 11.290 10.910 10.650 10.350 9.964 9.709 9.650 
0.035 + 0.001 12.020 11.620 11.350 11.040 10.654 10.379 10.344 



0.060 

0.024 
0.005 



" Sample stars used and derived fundamental physical parameters via IRFM. Apparent bolometric magnitudes (niBoi) have been computed according to Casagrande et al. 
(2006), where the absolute bolometric magnitude of the Sun Mboi,o = 4.74. For each star msoi is obtained using its bolometric flux and effective temperature and therefore it is 
already corrected for reddening, if present. Notice however that the observed magnitudes given here are not: before computing bolometric corrections, the observed magnitudes 
should be corrected using the corresponding E(B - V) given here. Errors have been computed as described in the text, without accounting for the uncertainty in E(B — V): 
changing it by ±0.01 would aff"ect T^s by approximately +50 K. 



